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PREFACE 


operations,  including  amphibious  assault  operations  for  Logistics-Over- 
the- Shore  (LOTS)  evolutions,  relies  extensively  on  the  utilization  of 
U.S.  Flag  commercial  shipping.  Since  the  mid  1960's  commercial  shipping 
has  been  steadily  shifting  towards  containerships,  Roll-on/Roll-off  (R0- 
RO)  ships  and  barge  ships  (e.g.,  LASH,  SEABEE) .  By  1985  as  much  as  85% 
of  U.S.  Flag  sealift  capacity  may  be  in  container  capable  ships  --  mainly 
non-selfsustaining  (NSS)  containerships.  Such  ships  cannot  operate  with¬ 
out  extensive  port  facilities.  Amphibious  assault  and/or  LOTS  operations 
are  usually  conducted  over  undeveloped  beaches  and  expeditious  response 
times  preclude  conventional  port  development.  Handling  of  containers  in 
this  environment  presents  a  serious  problem.  The  problem,  as  defined 
above,  is  addressed  in  the  overall  DOD  Over-the-Shore  Discharge  of  Cargo 
(OSDOC)  efforts  involving  developments  by  the  Army,  Navy  and  Marine  Corps 
Guiding  policy  is  documented  in  the  "DOD  Project  Master  Plan  for  Surface 
Container  Supported  Distribuiton  System"  and  the  OASD  I&L  system  defini¬ 
tion  paper  "Over-the-Shore  Discharge  of  Cargo  (OSDOC)  System." 


In  response  to  the  DOD  Master  Plan,  Navy  Operational  Requirement 
(0R-YSL03)  has  been  prepared  for  an  integrated  Container  Offloading  and 
Transfer  System  (COTS)  for  discharging  container  capable  ships  in  the 
absence  of  port  facilities.  The  COTS  Navy  Development  Concept  Paper 
(NDCP)  No.  YSL03  was  promulgated  July  1975  and  the  Navy  Material  Command 
tasked  with  development.  The  Naval  Facilities  Engineering  Command  has 
been  assigned  as  Principal  Development  Activity  (PDA)  with  the  Naval  Sea 
Systems  Command  assisting. 
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The  COTS  advanced  development  program  includes  the  ship  unloading 
subsystem,  the  ship-to-shore  subsystem  and  common  system  elements.  The 
ship  unloading  subsystem  includes:  (a)  the  development  of  Temporary  Con¬ 
tainer  Discharge  Facilities  (TCDF)  employing  merchant  ships  and/or  barges 
with  add-on  cranes  and  support  equipment  to  offload  non-selfsustaining 
(NSS)  containerships  alongside;  (b)  the  development  of  Crane  on  Deck  (COD) 
techniques  and  equipment  for  direct  placement  of  cranes  on  the  decks  of 
NSS  containerships  to  render  them  self sustaining  in  an  expedient  manner; 

(c)  the  development  of  equipment  and  techniques  to  offload  RO-RO  ships 
offshore  and  (d)  the  development  of  interface  equipment  and  techniques  to 
enable  ship  discharge  by  helicopters  (either  existing  or  projected  in 
other  development  programs).  The  ship-to-shore  subsystem  includes  the 
development  of  elevated  causeways  to  allow  cargo  handling  over  the  surfline 
and  development  of  self-propelled  causeways  to  transport  cargo  from  ships 
to  the  shores ide  interface.  The  commonality  subsystem  includes:  (a)  the 
development  of  wave  attenuating  Tethered  Float  Breakwaters  (TFB)  to  provide 
protection  to  COTS  operating  elements;  (b)  the  development  of  special 
cranes  and/or  crane  systems  to  compensate  for  container  motion  experienced 
during  afloat  handling;  (c)  the  development  of  transportability  interface 
items  to  enable  essential  outsize  COTS  equipment  transport  on  merchant 
ships  —  particularly  bargeships,  and  (d)  the  development  of  system  inte¬ 
gration  components  such  as  moorings,  fendering,  communications  and  services. 

This  report  addresses  the  progress  and  accomplishments  associated 
with  the  development  of  a  methodology  for  the  modeling  and  analysis  of  the 
complex  interactions  among  the  elements  of  the  ship  unloading  subsystem. 
Finally,  a  limited  sample  of  actual  test  data  is  used  to  validate  the  "COTS" 
computer  program  developed  for  this  study. 
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ABSTRACT 


A  strong  trend  in  the  commercial  shipping  industry  is  toward  con¬ 
tainerization  of  cargo,  a  system  that  utilizes  specialized  ships  and 
extensive  port  facilities.  With  dwindling  numbers  of  auxiliary  ships, 
the  military  is  forced  to  place  increased  reliance  upon  these  commer¬ 
cial  shipping  assets.  Since  military  operations  are  often  conducted 
over  unimproved  beaches,  an  integrated  container  handling  system  is 
required  to  support  a  broad  spectrum  of  supply  operations  in  the 
absence  of  improved  port  facilities. 

The  Navy's  Container  Offloading  and  Transfer  System  (COTS)  is 
under  development  to  meet  this  need.  Cranes  mounted  on  floating  plat¬ 
forms  will  be  used  to  lift  the  containers  from  the  cells  of  ships  and 
place  them  into  lighterage  for  transport  ashore.  When  this  operation 
is  conducted  in  the  presence  of  significant  wave  activity,  a  combina¬ 
tion  of  breakwaters  to  shelter  the  floating  vessels  and/or  motion  com¬ 
pensating  capability  of  the  cranes  may  be  required. 

A  methodology  is  developed  herein  for  the  evaluation  of  these  con¬ 
tainer  handling  systems.  The  wave  activity  is  treated  as  a  random 
forcing  function.  The  dynamics  and  interactions  of  the  system  elements 
are  described  by  appropriate  differential  equations  and  the  system 
measures  of  effectiveness  are  described  in  terms  of  the  responses  of 
the  dynamic  elements  and  the  overall  system  configuration.  The  computer 
model  is  compared  with  actual  test  results. 
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SECTION  I 


INTRODUCTION 


1 .1  BACKGROUND 

The  purpose  of  this  study  is  to  provide  a  methodology  for  the 
evaluation  of  alternative  container  offloading  system  configurations  in 
support  of  a  variety  of  military  amphibious  assault  operations.  This 
system  includes  the  interface  between  the  ships  involved,  their  dynam¬ 
ics,  the  offloading  dynamics,  and  the  dynamics  of  other  control  devices 
used. 

The  success  of  any  sustained  military  amphibious  assault  operation 
is  highly  dependent  upon  the  effectiveness  of  the  supporting  logistics 
resupply  activities.  In  the  event  of  a  military  contingency  involving 
the  landing  of  a  Marine  Amphibious  Force  (MAF)  of  approximately  50,000  men 
on  an  unimproved  beach,  the  resupply  pipeline  will  need  to  be  in  full 
operation  within  the  first  two  weeks  (approximately  two  weeks  of  sup¬ 
plies  go  ashore  with  the  MAF).^  Department  of  Defense  (DOD)  level 
planning  for  the  logistic  support  necessary  to  sustain  contingency  sea¬ 
lift  operations  relies  heavily  upon  the  utilization  of  U.S.  Flag  com¬ 
mercial  shipping  assets  to  augment  the  dwindling  number  of  Navy  and 
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Military  Sealift  Command  (MSC)  breakbulk  ships.  '  In  the  last  decade, 
the  dominant  trend  in  commercial  shipping  has  been  toward  high  speed, 
non-self-sustaining  container  ships,  and  by  1985,  as  much  as  85£  of  the 
U.S.  Flag  sealift  capabilities  may  be  in  these  container  capable 
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The  containership,  however,  requires  extensive  well -developed  port 

facilities  including  deep  draft  berths  and  large  industrial  gantry 

cranes.  The  response  time  and  general  nature  of  amphibious  assault 

operations  preclude  the  construction  of  such  port  facilities  in  their 

(41 

support  as  they  are  usually  conducted  over  undeveloped  beaches.  An 
integrated  and  balanced  cargo  handling  system  is  required  for  discharg¬ 
ing  container-capable  ships  in  support  of  a  broad  spectrum  of  military 
contingency  supply  operations  in  the  absence  of  improved  port  facil¬ 
ities.  The  system  could  also  find  peacetime  applications  for  emergency 
relief  operations  following  international  disasters  such  as  floods  or 
earthquakes. 

The  Naval  Facilities  Engineering  Command  (NAVFAC)  has  been  tasked 

as  the  Principal  Development  Activity  for  the  Container  Offloading  and 

(51 

Transfer  System  (COTS).  1  This  involves  the  development  of  the  re¬ 
quired  technology,  and  identification  and  procurement  of  the  hardware 
necessary  for  the  offloading  of  containers  from  ships  to  lighters 
(smaller  landing  craft  and  amphibious  or  air  cushion  vehicles)  in  an 
unprotected,  offshore  environment.  Current  development  efforts  are 
directed  at  the  application  of  some  form  of  crane  to  extract  the  con¬ 
tainers  from  the  cells  and  place  them  into  the  lighterage.  It  is  de¬ 
sirable,  from  the  standpoint  of  economy  and  availability,  to  utilize 
standard,  commercially  available,  mobile,  construction  cranes  for  this 
purpose.  However,  recent  tests  have  indicated  that  these  cranes,  de¬ 
signed  to  operate  on  land,  suffer  degraded  performance  when  employed  in 
unprotected,  offshore  environments. 
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This  environment  is  characterized  by  a  variety  of  natural  factors 
such  as  wind,  rain,  salt  spray,  and  wave  activity.  Wind  and  rain,  also 
experienced  by  shores ide  facilities,  do  not  pose  major  hurdles  except 
in  extreme  cases.  These  occur  with  such  small  likelihoods  that  oper¬ 
ations  can  be  terminated  temporarily  with  little  effect  on  the  overall 
system  productivity.  Much  work  has  been  done  in  the  protection  of 
machinery  from  the  harsh  marine  environment.  Significant  wave  activity, 
however,  occurs  with  such  regularity  in  offshore  environments  that  it 
cannot  be  dismissed  lightly.  The  random  wave  activity  at  sea  results 
in  irregular  motions  of  the  ships,  crane  platforms,  and  lighterage. 

These  wave  induced  motions  of  the  floating  crane  platform  are 
transmitted  through  the  crane  to  the  container.  Large  relative  dis¬ 
placements  between  the  container  and  the  adjacent  vessels  make  cargo 
transfer  both  difficult  and  dangerous.  Lateral  motions  near  the  natural 
frequency  of  the  suspended  load  are  amplified.  The  resulting  pendu- 
lation  may  introduce  high  side  loads,  and  thus  bending,  on  the  crane 
boom,  a  structural  member  designed  primarily  for  compression.  These 
effects  become  significantly  more  pronounced  as  the  sea  conditions 
worsen. 

A  variety  of  diverse  alternative  technologies  have  been  proposed 
to  ameliorate  the  adverse  effects  of  sea  induced  motions.  These  in¬ 
clude  the  use  of  larger  capacity  cranes  to  withstand  the  extra  loads, 
structural  modifications  or  add-on  features  to  existing  cranes  to 
improve  their  capabilities,  and  deployment  of  larger  quantities  of 
ships  and  cranes  to  offset  the  loss  of  productivity  resulting  from 


their  degraded  performance.  More  sophisticated  approaches  include  the 
development  of  automatic,  motion  compensating  cranes  and  readily  de¬ 
ployable  floating  breakwaters. 

The  motion  compensating  crane,  as  the  name  implies,  is  equipped 
with  a  means  of  sensing  the  sea  induced  motions  and  automatically 
raising  and  lowering  the  load  to  minimize  their  adverse  effects.  Oil 
exploration,  particularly  in  the  rough  North  Sea,  has  brought  about 
recent  commercial  interest  in  the  development  of  this  technology.  A 

variety  of  techniques  have  been  proposed  including  below  the  hook  de- 
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vices,'  '  hydraulic  cylinder  cable  tensioners,'  '  and  servo-controlled 

winches  with  control  strategies  ranging  from  simple  acceleration  feed- 

back^10^  to  complex  adaptive  controllers  approximating  trajectories 

derived  from  application  of  the  optimal,  quadratic,  linear  regulator 

results.  ^  ^ 

Breakwaters,  on  the  other  hand,  attack  the  problem  at  the  source. 

Breakwaters  attenuate  the  wave  energy  before  it  reaches  the  cargo 

handling  elements.  Fixed  breakwaters  have  been  used  for  many  years  for 
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harbor  protection  and  beach  erosion  control.'  '  Floating  breakwaters 

have  been  built  for  use  in  deep  water  applications  such  as  protecting 

offshore  oil  terminals. (^,14)  Aether  operating  on  the  principle  of 

energy  reflection  or  dissipation,  they  tend  to  be  large,  massive 

structures.  Current  development  efforts  are  directed  at  improving  the 

effectiveness  of  floating  breakwaters  that  are  easily  transportable 
(31 

and  deployable.'  ' 

In  addition  to  the  above  technologies,  a  variety  of  basic  system 
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configurations  have  been  proposed'  ’  '  including  the  Crane  on  Deck  (COD) 


and  the  Temporary  Container  Discharge  Facility  (TCDF).  The  COD  ap¬ 
proach  suggests  that  mobile  cranes  can  be  placed  on  the  deck  of  the 
container  ships  to  give  them  self-off loading  capability.  The  TCDF  ap¬ 
proach  envisions  a  dedicated  platform,  either  a  ship  or  a  large  barge, 
that  can  be  equipped  with  the  necessary  cranes  and  personnel  to  remain 
in  the  objective  area  and  act  as  an  offshore  floating  pier. 

Because  of  the  diverse  nature  of  the  technologies  that  affect  the 
system  design  and  performance,  the  complex  nature  of  the  interactions 
between  the  various  system  elements,  and  the  large  number  of  possible 
system  configurations,  there  has  been  no  quantitative  means  of  compar¬ 
ing  and  evaluating  alternative  system  design  candidates. 

1.2  PROBLEM  DESCRIPTION 

The  problems  faced  by  designers  of  the  Container  Offloading  and 
Transfer  System  (COTS)  can  be  divided  into  two  major  areas.  The  first 
is  the  description  of  the  complex  interactions  between  the  system 
elements  and  how  these  interactions  affect  the  system  measures  of  ef¬ 
fectiveness.  The  second  is  the  determination  of  the  optimal  allocation 
of  available  resources  to  maximize  the  overall  system  performance. 

The  COTS  problem  is  not  unique,  but  a  member  of  a  larger,  more 
general  class  of  problems.  Many  complex  systems  of  dynamic  elements 
are  driven  by  random  inputs  where  system  performance  is  measured  in 
terms  of  descriptive  statistics  of  the  system's  response.  In  designing 
to  required  performance  criteria,  a  tradeoff  may  exist  between  re¬ 
duction  of  the  undesirable  aspects  of  the  random  input  or  modification 
of  the  dynamic  characteristics  of  the  system  itself. 


Other  examples  of  such  systems  include  high  speed  rail  transport¬ 
ation  systems  and  short  take-off  and  landing  (STOL)  air  terminal  de¬ 
sign.  In  the  case  of  rail  transportation,  the  tradeoff  might  be  bet¬ 
ween  the  costs  of  straightening  the  rails  to  various  tolerances  versus 
implementation  of  anticipatory  control  suspension  hardware  on  the 
trains  themselves.  Performance  criteria  might  be  mean  square  linear 
acceleration  or  noise  level  experienced  by  the  passengers.  In  the  case 

of  STOL  air  terminals,  the  costs  of  wind  screens  may  be  traded  off 
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against  more  sophisticated  control  strategies  for  the  aircraft  ' 
where  expected  touchdown  accuracy  might  be  the  index  of  performance. 

The  COTS  design  problem  is  similarly  organized.  The  random  input 
from  the  sea  can  be  reduced  by  the  addition  of  breakwaters.  On  the 
other  hand,  a  variety  of  system  configurations  and  crane  designs  that 
ameliorate  the  adverse  effects  of  this  input  are  available.  Overall 
cargo  handling  productivity  and  system  cost  are  governed  by  the  design 
of  the  system  as  a  whole. 

1.3  OBJECTIVE 

The  objective  of  this  effort  is  to  develop  a  methodology  by  which 
alternative  Container  Offloading  and  Transfer  Systems  (COTS)  can  be 
compared  and  evaluated  with  respect  to  their  productivity  and  cost. 

This  will  include  the  development  of  an  analytical  model  to  predict 
the  system  cost  and  performance  in  terms  of  descriptive  parameters  of 
the  system  elements.  With  this  model,  optimization  studies  will  be 
conducted  on  certain  parameters  of  interest.  Sensitivity  analyses  can 
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be  performed  on  the  optimal  solutions  with  respect  to  system  parameters 
and  constraints  to  help  identify  profitable  areas  for  further  research 
and  development. 

1.4  SYSTEM  ELEMENTS 

The  Container  Offloading  and  Transfer  System  (COTS)  is  composed  of 
eight  basic  elements  -  sea  state,  breakwater,  crane  platform,  container 
ship,  lighterage,  crane,  cargo,  and  personnel.  These  are  described 
briefly  in  Table  1.1  along  with  typical  parameters  commonly  used  in 
their  quantification.  These  elements  are  also  arranged  in  an  inter¬ 
action  matrix  (Figure  1.1),  providing  a  visual  representation  of  the 
kinds  of  relationships  that  must  be  identified  in  order  to  effectively 
model  the  system. 


Sea  Stat 

Crane  PI 
Containe 

Lighters 

Crane 

Cargo 


Figure  1.1.  COTS  Interaction  Matrix 
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TABLE  1.1  DESCRIPTION  OF  SYSTEM  ELEMENTS 


SYSTEM  ELEMENTS  BRIEF  DESCRIPTION  TYPICAL  PARAMETERS 


Sea  State 

Measure  of  wave  activity 

Significant  wave  height 
power  density  spectrum 

Breakwater 

Wave  attenuating  device 

Transmission  co¬ 
efficients,  length 

Crane  Platform 

Floating  vessel  on 
which  the  crane  is  placed 

Length,  breadth, 
draft,  weight,  hull 
design 

Container  Ship 

Large  oceangoing  vessels 
for  transporting  contain¬ 
erized  cargo 

Length,  breadth, 
draft,  weight,  hull 
design 

Lighterage 

Small  floating  vessels 
used  for  transporting 
personnel  and  logistics 
from  ships  to  shore 

Length,  breadth, 
draft,  weight,  hull 
design 

Crane 

Mechanical  structure 

for  lifting  and  moving 
heavy  loads 

Lift  capacity,  reach, 
horsepower,  stiffness, 
design 

Cargo 

Supplies  to  be  handled 
by  the  system 

Weight,  dimensions, 
impact  resistance 

Personnel 

Collection  of  Individuals 
necessary  to  operate  the 
system 

Training,  motivation 

With  the  system  bounded  to  include  these  basic  elements,  it  is 
worthwhile  to  consider  the  environment  in  which  the  system  will  operate 
and  the  influence  it  may  exert  on  the  system  elements.  Such  factors  as 
weather  conditions  other  than  wave  activity,  personnel  training  and 
motivation,  availability  of  lighterage,  and  the  like  must  be  considered. 
These  may  vary  greatly  from  one  scenario  to  another,  but  for  the  pur¬ 
pose  of  this  study,  nominal  conditions,  based  on  past  Department  of 
Defense  (DOD)  experience,  will  be  assumed. ^ 

1.5  MEASURES  OF  EFFECTIVENESS 

i 

There i are  a  number  of  system  measures  of  effectiveness  that  re¬ 
flect  DOD  goals  and  objectives.  They  include  cost,  productivity, 
transportability,  reliability,  response  time,  personnel  requirements, 
and  safety.  Of  these,  cost  and  productivity  are  the  more  easily  quan¬ 
tifiable.  Further,  the  others  can  be  represented  in  terms  of  these 
two.  For  example,  a  Temporary  Container  Discharge  Facility  (TCDF), 

consisting  of  a  crane  mounted  on  a  barge,  can  be  towed  at  a  maximum 
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speed  of  five  knots.  '  It  could,  however,  be  loaded  onto  a  SEABEE 
ship  (a  commercial  vessel  designed  to  carry  larae  barges)  and  trans¬ 
ported  at  20  knots  to  the  objective  area.^3,16^  Because  there  is  a 
limited  number  of  these  ships,  space  on  them  would  be  worth  a  premium 
during  a  military  contingency,  but  a  dollar  value  can  be  placed  on  its 
use.  This  would  provide  a  basis  to  assign  a  cost  to  this  option  for 
evaluation  purposes.  Based  upon  this  line  of  reasoning,  cost  and  pro¬ 
ductivity  are  chosen  as  useful  parameters  for  the  purpose  of  alter¬ 
native  evaluations. 


1.6  SYSTEM  SYNTHESIS 


In  addition  to  the  descriptive  parameters  of  the  system  elements, 
system  configuration  information  is  needed  to  describe  the  dynamic  in¬ 
teractions  of  the  system  elements.  System  configuration  includes  in¬ 
formation  relating  to  system  operation,  arrangement  of  the  elements, 
the  mooring  plan,  and  the  like. 

In  order  to  insure  the  broad  capability  of  a  model  to  evaluate  a 
variety  of  COTS  alternatives,  other  miscellaneous  inputs  may  be  re¬ 
quired  such  as  fixed  cost  information  or  baseline  productivity  data 
that  are  particular  to  the  configurations  under  investigation. 

The  basic  structure  of  the  model  is  shown  in  Figure  1.2.  The 
model  is  divided  into  two  major  areas  -  system  dynamics  and  system 
evaluation.  The  system  dynamics  model  addresses  the  system  inter¬ 
actions  that  can  be  described  by  the  appropriate  application  of  phys¬ 
ical  laws  of  nature.  This  includes  the  calculation  of  transmitted 
wave  activity,  ship  and  platform  motions,  crane  structural  require¬ 
ments,  and  control  system  response.  The  system  evaluation  model,  how¬ 
ever,  is  less  well  defined.  It  includes  the  estimation  of  component 
costs  and  system  productivity  based  upon  outputs  from  the  system  dy¬ 
namics  model.  The  relationships  necessary  for  system  evaluation  will 
be  developed  empirically,  extrapolated  from  existing  data,  and 
structured  in  a  manner  that  appears  most  reasonable  to  experts  familiar 
with  the  Container  Offloading  and  Transfer  System  (COTS). 

A  system  synthesis  flow  chart  (Figure  1.3)  shows  the  interactions 
between  the  various  system  elements  and  the  contributions  to  the  system 
measures  of  effectiveness.  It  provides  a  visual  representation  of  the 


SYSTEM  CONFIGURATION 


Figure  1.2.  Basic  Structure  of  the  COTS  Model 


Figure  1.3.  System  Synthesis  Flow  Chart 
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relationship  between  the  system  dynamics  model  and  the  system  eval¬ 
uation  model.  The  ancillary  interactions  between  personnel,  lighter¬ 
age,  and  cranes  are  also  shown,  although  they  will  be  implicitly  in¬ 
cluded  in  the  empirical  relationships  for  productivity. 


SECTION  II 


METHODOLOGY 


2.1  APPROACH 

In  developing  the  models  necessary  to  effectively  describe  the 
system  interactions,  a  number  of  simplifying  assumptions  are  made. 

This  is  done  in  accordance  with  conventional  engineering  practice  in 
order  to  enhance  the  mathematical  tractability  of  the  problem  and  to 
make  use  of  a  wide  body  of  powerful  analytical  techniques.  The  advan¬ 
tages,  consequences,  and  limitations  of  these  assumptions  are  addressed. 

The  system  dynamics  are  described  by  differential  equations  that 
are  assumed  to  be  linear  or  linearizable  about  an  operating  point. 

The  coefficients  of  these  equations  are  constant.  Although  few,  if  any, 
physically  realizable  systems  can  be  exactly  described  by  linear  dif¬ 
ferential  equations  with  constant  coefficients,  many  can  be  approximated 
in  this  manner  over  limited  operating  ranges. 

An  important  characteristic  of  linear  systems  is  frequency  pre¬ 
servation.  A  sinusoidal  input  of  a  given  frequency  results  in  a  sin¬ 
usoidal  output  at  the  same  frequency,  in  which  only  amplitude  and  phase 
are  modified. For  this  reason,  a  useful  characterization  of 
linear  systems  is  by  their  frequency  response  operators  or  transfer 
functions.  These  transfer  functions  can  be  obtained  by  Laplace  trans¬ 
formation  of  the  governing  differential  equations  or  by  experimental 
methods  on  the  hardware  itself. 
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The  governing  differential  equations  that  describe  the  motion  of 

the  sea  surface  can  also  be  simplified.  The  fluid  is  assumed  to  be 

inviscid,  the  Bernoulli  equation  is  linearized,  the  boundary  conditions 

of  the  free  surface  are  satisfied  at  its  mean  value  and  not  along  the 

wave  profile,  and  the  pressure  at  the  surface  is  assumed  to  be  con- 
(19) 

stant.  '  This  results  in  a  sinusoidal  motion  of  the  sea  surface. 

The  variation  in  the  level  of  the  sea  surface  at  any  point  results 
from  the  summation  of  contributions  from  many  relatively  unrelated 
disturbances.  This  collection  of  sinusoids  of  varying  amplitudes  and 
phase  angles  is  assumed  to  be  an  ergodic,  Gaussian  random  process. 
Statistical  analysis  of  actual  wave  records  indicates  that,  after 
elimination  of  the  linear  trends,  the  Gaussian  assumption  is  quite 
reasonable/20^ 

This  linearized,  Gaussian  model  of  a  seaway  is  limited  by  the 
fact  that  it  is  not  valid  for  high  or  extreme  sea  conditions.  The 
linearization  introduces  errors  into  the  profile  of  the  higher  waves. 

In  a  strictly  Gaussian  process,  extremely  high  waves,  although  im¬ 
probable,  would  be  possible,  but  waves  are  limited  in  height  by 
breaking. 

For  the  purpose  of  the  COTS  modeling,  these  do  not  represent 
objectionable  limitations.  The  operational  requirements  of  moderate 
sea  conditions  are  well  within  the  capabilities  of  the  linearized 
models.  Under  extreme  sea  or  storm  conditions  where  this  model  is  less 
reliable,  the  container  handling  operations  would  not  be  conducted. 


16 


The  advantage  of  the  linearized  Gaussian  model  is  that  a  wide  body 
of  powerful  mathematical  tools  are  available  for  the  analysis  of  systems 
with  these  forcing  functions.  A  zero  mean  random  process  of  this  type 
can  be  completely  described  by  its  power  density  spectrum.  This  is  a 
measure  of  the  energy  density  of  the  random  process  as  a  function  of 
frequency.  This  is  obtained  by  taking  the  Fourier  transform  of  the  auto 
correlation  function  or  may  be  determined  by  experimentally  measuring, 
filtering,  and  squaring  a  time  history  of  the  process.  A  narrowband 
random  process  is  one  with  a  peaked  power  density  spectrum,  indicating 
that  the  majority  of  its  energy  is  concentrated  in  a  narrow  frequency 
range. 

It  has  been  shown  that  a  linear  system,  subjected  to  a  Gaussian 

random  input,  has  a  Gaussian  output. Lord  Rayleigh  demonstrated 

that  the  peaks  of  a  narrow  band,  Gaussian  process  are  distributed 

(17  22) 

according  to  a  Rayleigh  probability  density  function.'  ’  '  The  power 

density  spectrum  of  the  output,  So(u),  of  a  linear  system  can  be  calcu¬ 
lated  from  the  power  density  spectrum  of  the  input,  S. («),  and  the  sys¬ 
tem's  complex  transfer  function,  H(jw),  by  the  relationship: 

SQ(w)  =  H( jw)  H*( jco)  SiU),(17’21) 

where  the  asterisk  indicates  the  complex  conjugate. 

The  standard  deviation,  and  therefore  a  complete  description  of 
the  Gaussian  distribution  of  a  zero  mean  random  process,  can  be  calcu¬ 
lated  from  its  power  density  spectrum. 
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where  a  is  the  standard  deviation  and  S(w)  is  the  single-sided  power 
density  spectrum. 

The  system  dynamics  model  is  based  on  the  above  analytical  tech¬ 
niques.  The  various  system  elements  or  groups  of  elements  are  described 
by  their  power  density  spectra.  The  standard  deviations  of  some  of  the 
outputs  will  be  used  by  the  system  evaluation  model  in  order  to  calcul¬ 
ate  expected  values  of  those  outputs  that  affect  the  cost  and  produc¬ 
tivity  of  the  system. 

2.2  SYSTEM  PARAMETERS 

In  light  of  the  above  assumptions,  it  is  useful  to  examine  the 
system  elements  and  the  types  of  parameters  that  are  descriptive  of 
their  role  in  the  system. 

2.2.1  Sea  State  Parameters 

The  sea  state  has  been  discussed  above.  It  will  be  considered  as 
a  narrowband,  Gaussian,  random  process.  Sea  state  will  be  described  by 
its  power  density  spectrum.  The  sea  is  a  continuously  changing  source 
of  activity  driven  by  wind  velocity,  duration,  fetch,  bottom  contour, 
and  the  like.  Because  of  the  complex  nature  of  the  wave  generation  pro¬ 
cess,  the  power  density  spectrum  changes  with  time.  Oceanographers 
have  examined  the  growth  of  wave  spectra  as  the  seas  develop  under  a 
variety  of  conditions.  For  a  given  set  of  wind  conditions,  a  steady 
state  wave  spectrum  is  reached,  beyond  which  no  further  increase  in 
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energy  is  observed.  This  is  called  a  fully  developed  sea.  Mariners, 
on  the  other  hand,  have  a  relatively  simple,  although  operationally 
effective  method  of  describing  the  wave  activity  associated  with  sea 
states.  This  is  significant  wave  height,  defined  as  the  average  of  the 
one  third  largest  waves  as  measured  from  their  peak  to  trough.  For  the 
purposes  of  the  COTS  modeling,  an  analytical  form  of  the  power  density 
spectrum  of  a  fully  developed  sea  that  can  be  defined  in  terms  of  sig¬ 
nificant  wave  heights  will  be  used. 

2.2.2  Breakwater  Parameters 

There  are  two  basic  types  of  breakwater  -  fixed  and  floating. 

Most  fixed  breakwaters  provide  a  rigid,  full  depth  barrier  to  the  pro¬ 
gress  of  the  waves.  Floating  breakwaters  provide  a  movable,  partial 
barrier  and  operate  by  a  combination  of  energy  reflection  and  dissipa¬ 
tion.  Figure  2.1  shows  several  breakwater  configurations. 

With  the  great  variety  in  the  designs  of  breakwaters,  parameters 
descriptive  of  the  physical  characteristics  would  be  difficult  to 
handle  effectively  in  the  model.  Since  the  primary  interaction  of  the 
breakwater  is  with  the  sea  state,  parameters  associated  with  the  trans¬ 
mitted  wave  height  would  be  most  useful.  A  casual  examination  of  typ¬ 
ical  floating  breakwater  wave  transmission  coefficients  suggests  that 
breakwaters  have  frequency  characteristics  similar  to  those  of  low  pass 
filters.  The  longer  wave  length,  lower  frequency  waves  are  transmitted 
with  very  little  attenuation.  The  shorter,  higher  frequency  waves,  how¬ 
ever,  are  significantly  reduced.  This  motivates  a  description  of  break¬ 
water  transmission  coefficients  as  a  function  of  cut  off  frequency  and 
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Figure  2.1.  Typical  Breakwater  Configurations 
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high  frequency  attenuation. 

The  Tethered  Float  Breakwater  (TFB)  is  currently  under  development 

at  the  Naval  Ocean  Systems  Center  and  Scripps  Institution  of  Ocean- 
(23) 

ography  for  possible  COTS  applications.  A  comparison  of  trans¬ 
mission  coefficients  and  cost  for  several  alternative  TFB  designs  in¬ 
dicate  that  cost  may  also  be  approximated  by  a  function  of  the  para¬ 
meters  that  define  the  transmission  coefficient. 

2.2.3  Ship,  Lighterage,  and  Platform  Parameters 

The  Container  Offloading  and  Transfer  System  includes  three  dif¬ 
ferent  types  of  floating  vessels  -  the  crane  platform,  the  container- 
ship,  and  the  lighterage.  All  of  these  contribute  to  the  cargo  han¬ 
dling  problem  by  virtue  of  their  response  to  the  sea  state.  The  mo¬ 
tions  of  the  crane  platform  result  in  amplified  motions  of  the  sus¬ 
pended  cargo,  inducing  dynamic  loads  in  the  crane  structure.  Relative 
motions  between  the  other  ships  and  the  cargo  result  in  degraded  cargo 
handling  productivity. 

A  floating  vessel,  assumed  to  be  rigid,  has  six  degrees  of 
freedom  -  heave,  surge,  sway,  roll,  pitch,  and  yaw.  The  ratio  of  the 
amplitude  of  the  response  of  a  ship  to  the  amplitude  of  the  incident 
waves  as  a  function  of  frequency  is  called  the  response  amplitude  oper¬ 
ator  (RAO).  It  has  been  shown  that  these  RAO's  are  linear  operators 
for  wall-sided  ships  oscillating  in  light  to  moderate  seas.^24,25,26^ 
The  RAO's  for  the  six  degrees  of  freedom  of  a  ship  may  be  calculated 

experimentally  by  recording  the  response  of  a  scale  model  in  a  wave 
(27) 

tank.  Analytical  techniques  have  been  developed  to  calculate  the 
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RAO's  through  solution  of  linearized  differential  equations.'  *  ' 

The  information  necessary  to  describe  the  ship  includes  length,  breadth, 
draft,  hull  design,  weight,  moments  of  inertia,  natural  roll  period, 
water  depth,  and  direction  of  the  incident  waves. 

2.2.4  Crane  Parameters 
2. 2. 4.1  Conventional  Cranes 

A  crane  is  a  mechanical  device  designed  to  hoist  and  move  a  load 
in  a  controlled  manner.  Cranes  come  in  a  variety  of  shapes  and  sizes 
(Figures  2.2).  They  can  be  divided  into  two  major  types  -  those  with 
traveling  beams  such  as  overhead  bridge  and  gantry  cranes  and  those 
with  revolving  jibs  such  as  mobile  construction  cranes  and  shipboard 
cargo  handling  cranes.  Gantry  cranes  have  been  employed  effectively 
on-board  ships  for  the  offloading  of  containerized  cargo.  These  cranes, 
however,  are  specially  designed  for  each  ship  and  crane  rails  must  be 
provided  on  the  deck  of  the  ship.  Only  a  limited  number  of  container 
ships  are  equipped  with  these  features'  '  and  they  may  not  be  readily 
available  on  short  notice.  The  response  of  a  COTS  to  a  military  re¬ 
quirement  precludes  the  use  of  such  cranes.  Mobile  construction 
cranes  of  the  revolving  boom  type  have  sufficient  lift  and  reach  capa¬ 
bilities  to  meet  the  COTS  requirements  and  are  readily  available  from 
the  commercial  construction  industry.  ’  '  It  has  been  determined 

that  these  cranes  represent  viable  candidates  for  the  COTS  crane 
application.  It  is  desirable  to  analyze  their  capabilities  in  the  com¬ 
mercially  available  state,  with  minor  modifications  such  as  power 
taglines  or  a  Rider  Block  Tagline  System  (RBTS),  and  with  major 
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Figure  2. 2. a.  Traveling  Beam  Cranes 


modifications  such  as  automatic  vertical  motion  compensation. 

The  primary  parameters  that  describe  COTS  cranes  are  their  type, 
lift  capacity,  and  reach  capability.  Mobile  construction  cranes  are 
basically  of  two  major  types  -  truck  mounted  and  crawler  mounted, 
although  special  modifications  like  the  Ringer^32^  and  the  Skyhorse^33^ 
are  available  on  certain  cranes  to  greatly  increase  their  lift  -  reach 
capability. 

The  position  of  a  load  being  handled  by  a  revolving  jib  crane  is 
controlled  by  a)  raising  and  lowering  the  hook,  called  hoisting, 
b)  changing  the  vertical  angle  of  the  boom,  called  luffing,  and  c)  ro¬ 
tating  the  crane  structure  about  its  vertical  axis,  called  slewing. 

Lift  capacity  is  the  load  that  the  crane  can  safely  handle.  Reach  is 
the  radius  at  which  it  can  handle  the  load. 

Typical  lift  -  reach  characteristics  of  a  commercially  available 
crawler  crane  are  displayed  in  Figure  2.3.  It  can  be  seen  that  the 
lifting  capacity  is  significantly  reduced  as  the  operating  radius  is 
increased,  where  overturning  of  the  crane  becomes  the  primary  limiting 
factor. 

Other  characteristics  of  a  crane  that  influence  the  system  perfor¬ 
mance  are  the  crane's  operating  speeds,  the  vertical  stiffness  of  the 
crane  structure,  and  its  pendulation  characteristics.  The  maximum 
crane  operating  speeds  place  a  practical  limit  on  the  maximum  container 
offloading  rate.  The  combination  of  the  vertical  stiffness  of  the 
crane  structure  and  the  mass  of  the  lifted  load  constitutes  a  dynamic 
system  that  may  amplify  the  effective  load  on  the  crane  at  certain 


frequencies  and  contribute  to  the  maximum  impact  loads  experienced  dur¬ 
ing  lift  off  of  the  containers  from  the  deck.^’3^’3^  The  natural 
period  of  the  pendulus  load  suspended  from  the  boom  tip  is  a  function  of 
the  length  of  the  pendulum.  This  lightly  damped  system  can  be  excited 
by  the  motions  of  the  crane  platform.  The  resulting  pendulation  if  un¬ 
restrained,  may  introduce  excessive  sideloads  on  the  crane  boom,  a 
structural  member  designed  primarily  for  axiul  compression.  Development 
efforts  have  been  directed  at  reducing  this  effect.  Hydraulically  pow¬ 
ered  taglines  reeved  from  winches  near  the  boom  base  to  the  load  have 
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been  employed  to  introduce  damping,  '  and  the  Rider  Block  Tagline  Sys¬ 
tem  (RBTS)  has  been  developed  to  reduce  the  adverse  effects  of  pendu¬ 
lation. (38) 


2. 2. 4. 2  Rider  Block  Tagline  System 

The  Rider  Block  Tagline  System,  shown  in  Figure  2.4,  consists  of  a 
rider  block,  outriggers,  tagline,  and  a  rider  block  line  supplementing 
the  conventional,  commercial  crane's  hardware.  The  load  lines  pass  over 
sheaves  in  the  rider  block.  The  horizontal  and  vertical  positions  of 
the  rider  block  are  controlled  by  the  operator  with  the  taglines  and  the 
rider  block  line,  respectively.  This  allows  the  operator  to  control  the 
pendulation  length  and  load  position  simultaneously.  Figure  2.5  shows 
the  rider  block  in  greater  detail.  It  is  fully  articulated  so  that  it 
can  conform  to  the  line  of  action  of  each  wire. 

Actuation  of  the  taglines  or  the  rider  block  line  results  in  a  non¬ 
linear  response  of  the  load  that  may  require  a  degree  of  operator/signal 
man  proficiency  to  insure  a  smooth  operation.  Figure  2.6  depicts  the 
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Figure  2.5.  Rider  Block  Components 
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nature  of  this  response  for  equal  increments  of  tagline  and  rider  block 
line  position. 

The  principal  advantage  of  the  RBTS  is  that  it  redirects  the  forces 
caused  by  pendulation  and  ship  roll  from  sideloading  of  the  boom  into 
differential  tagline  tensions.  The  RBTS  also  reduces  the  effective 
length  of  the  pendulum,  providing  a  more  controllable  load  handling 
operation. 

2. 2. 4. 3  Motion  Compensating  Cranes 

Motion  compensating  cranes  represent  an  extension  of  conventional 
crane  technology  by  incorporating  motion  sensing  and  automatic  control 
circuitry  with  the  hoisting  machinery  to  reduce  the  relative  motion  bet¬ 
ween  the  cargo  and  the  lighterage  or  containership.  Development  of  an 
effective  motion  compensating  crane  include  sensing  hardware,  a  control 
strategy,  and  a  means  of  implementing  the  control. 

A  number  of  reports  have  addressed  the  problem  of  sensing  the 
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relative  motion  between  ships  at  sea.'  ’  '  Solutions  include: 

a)  Accelerometers  mounted  in  the  load  block,  crane  boom 
tip,  and  on  the  deck  of  the  adjacent  ship. 

b)  Laser  or  radar  units  mounted  on  the  container  spreader 
bar  directed  downward. 

c)  A  system  of  triangulating,  angle  sensing,  taut  wires 
from  the  crane  platform  to  the  load  and  the  adjacent 
ship. 

d)  A  cable  reel  mounted  on  the  container  spreader  with  a 
line  that  can  be  hand  held  by  a  cargo  handler  on  or 
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temporarily  attached  to  the  adjacent  craft  to  provide 

a  direct  measure  of  relative  vertical  displacement. 

The  first  control  strategies  examined  in  detail  were  those  that 

synchronized  the  load  with  the  deck  of  the  ship  on  which  it  was  to  be 
( 41 *42 1 

placed.'  ’  They  included  variations  of  acceleration,  velocity,  and 
position  feedback  systems,  but  all  involved  consumption  of  large  amounts 
of  horsepower,  resulting  in  massive  hoist  machinery  designs.  Like  most 
servo  systems,  the  bandwidth  of  the  response  is  reduced  as  the  masses  of 
the  moving  parts  increase.  The  necessary  horsepower  could  be  obtained 
from  diesel  engines,  but  applying  it  at  the  winch  drum  in  a  controlled 
manner  is  a  more  formidable  problem. 

(41  431 

Several  methods  have  been  proposed  including:'  *  ' 

a)  Variable  pressure  clutch  and  brake  systems. 

b)  Electric  motor  and  generator  systems  of  the 
Ward  Leonard  type. 

c)  Hydraulic  motors  driven  by  variable  displacement 
pumps . 

The  clutch  -  brake  system  Is  responsive  but  has  severe  wear  problems 
over  extended  periods  of  time.  The  electric  motor  -  generator  systems 
have  time  constants  that  are  highly  dependent  on  the  effective  inertia 
of  the  system,  and  for  the  COTS  loads,  are  not  sufficiently  responsive. 
The  hydraulic  systems  are  commercially  available,  are  easily  adaptable 
to  diesel  drives,  and  have  adequate  response  characteristics  that  are 
relatively  insensitive  to  load. 

Adaptive  control  strategies  have  been  examined  to  reduce  the  horse- 
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power  requirements.^^  The  objective  is  to  land  the  cargo  near  the 
crest  of  a  wave  where  the  relative  velocities  will  be  small,  and  there¬ 
fore,  the  synchronization  energy  minimized.  This  work,  however,  is  still 
in  the  theoretical  stages,  and  practical  implementation  of  the  strategy 
in  the  presence  of  structural  resonances  and  human  operator  interfaces, 
has  yet  to  be  demonstrated. 

An  alternative  motion  compensating  strategy  that  is  within  the 
range  of  current  commercial  capabilities  will  be  considered  here.  It  is 
the  platform  stabilized  motion  compensating  crane.  As  the  crane  plat¬ 
form  rolls  and  heaves,  the  load  experiences  random  vertical  displacements 
in  inertial  space.  A  control  system  based  on  velocity  feedback  could  be 
employed  to  hoist  or  lower  the  load  line  in  order  to  reduce  this  motion. 
This  technique  is  examined  in  greater  detail  in  Appendix  A. 2.  Such  a 
system  can  be  represented  by  transfer  functions  and  analyzed  and  stabi¬ 
lized  by  classical  control  theory  techniques.  The  transfer  function 
between  the  vertical  motion  of  the  crane  boom  and  the  load  has  the  char¬ 
acteristics  of  a  high  pass  filter.  That  is,  low  frequency  excursions  of 
the  crane  boom  can  be  compensated,  but  the  higher  frequency  motions 
would  be  transmitted  to  the  load  with  very  little  attenuation.  The  low 
frequency  attenuation  and  the  bandwidth  of  the  system  can  be  directly 
converted  into  horsepower  and  hardware  requirements  for  any  given  sea 
conditions.  This  can  be  reflected  in  the  system  cost,  and  the  re¬ 
duction  in  relative  motion  of  the  cargo  can  be  translated  into  increased 
productivity. 
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2.3  MODEL  SYNTHESIS 

With  all  of  the  relevant  dynamic  characteristics  of  the  elements 
described  by  transfer  functions  (Energy  Transmission  Ratio  (ETR), 
Response  Amplitude  Operator  (RAO),  and  the  like),  the  input  -  output 
characteristics  of  the  dynamic  model  can  be  determined.  The  dynamic 
model  is  a  single  input  -  multiple  output  system.  The  single  input  is 
the  wave  activity  or  sea  state  described  in  terms  of  its  power  density 
spectrum. 

By  the  appropriate  arrangement  of  the  transfer  functions  of  the 
dynamic  elements,  the  power  density  spectra  of  the  desired  outputs  can 
be  calculated.  These  include: 

a)  the  transmitted  sea  state, 

b)  the  motions  of  the  vessels, 

c)  the  dynamic  loadings  on  the  crane, 

d)  the  relative  displacement  between  the  cargo  and 
adjacent  ships,  and 

e)  the  hoist  line  velocities  of  the  motion  compensating 
machinery. 

These  are  all  spectra  of  Gaussian  random  processes,  whose  variance 

is  the  integral  from  zero  to  Infinity  of  its  spectrum.  The  square  root 

of  the  variance  is  the  standard  deviation.  It  can  be  shown  that  the 

significant  amplitude,  defined  as  the  average  of  the  one-third  highest 

excursions  of  the  random  process,  is  approximately  equal  to  twice  the 
(221 

standard  deviation.'  '  The  average  of  the  1 /1 000th  highest  excursions, 
often  used  as  a  measure  of  the  maximum  value,  equals  3.7  times  the 
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standard  deviation.^0,  The  significant  and  maximum  values  of  the 
system  responses  can  be  calculated  by  the  above  method  and  supplied  to 
the  system  evaluation  model. 

The  system  evaluation  model  will  calculate  a  cost  and  productivity 
for  the  total  system  using  the  nominal  cost  input  data,  breakwater  cost 
estimates  based  on  size  and  frequency  characteristics,  crane  cost  esti¬ 
mates  based  on  crane  dynamic  loadings  and  horsepower  requirements,  nom¬ 
inal  crane  characteristics,  and  significant  relative  displacements  bet¬ 
ween  the  load  and  the  adjacent  vessels. 

Figure  2.7  is  an  expanded  system  synthesis  flow  chart  showing  the 
block  diagram  of  the  system  dynamics  transfer  functions  and  the  flow  of 
information  to  the  system  evaluation  models. 

2.4  EVALUATION  TECHNIQUES 

The  COTS  model  can  be  used  in  a  variety  of  ways.  The  most  direct 
is  in  the  evaluation  of  alternative  COTS  configurations.  In  this  mode, 
the  alternative  systems  are  completely  described  by  the  user  and  the 
model  would  evaluate  the  cost  and  productivity  of  each  alternative. 

Another  use  might  be  to  generate  productivity  versus  sea  state 
curves  for  certain  alternative  designs.  This  information  could  be  used 
to  reevaluate  the  sea  state/productivity  requirements  currently  placed 
on  the  system  or  used  in  conjunction  with  sea  state  probability  infor¬ 
mation  to  predict  expected  container  offloading  capabilities  over  ex¬ 
tended  periods  of  time. 

Tradeoff  studies  between  competing  technologies  can  also  be  con¬ 
ducted  with  this  model.  The  primary  tradeoff  in  this  COTS  analysis  is 


between  motfon  compensating  cranes  and  breakwaters  as  a  means  of  im¬ 
proving  the  productivity  of  the  system.  The  breakwater  affects  the 
cargo  motion  as  a  low  pass  filter  while  the  motion  compensating  crane 
appears  as  a  high  pass  filter.  The  costs  of  both  breakwaters  and  motion 
compensating  cranes  increase  rapidly  as  their  bandwidths  are  broadened 
to  include  the  majority  of  the  sea  induced  cargo  and  ship  motions.  It 
is,  therefore,  reasonable  to  investigate  whether  or  not  there  is  some 
optimum,  intermediate  design  of  breakwaters  and  motion  compensating 
cranes  that  will  provide  a  required  level  of  productivity  at  minimum 
cost.  A  steepest  descent  algorithm  with  numerical  partial  deviations  of 
cost  and  productivity  with  respect  to  descriptive  parameters  of  the 
crane  and  breakwater  design  or  a  parameter  search  technique  could  be 
employed  in  the  solution  of  this  minimum  cost  tradeoff  problem. 

The  model  would  also  be  useful  in  conducting  a  sensitivity  analysis 
of  a  particular  configuration's  cost  and  productivity  to  all  of  the 
descriptive  parameters  of  the  system  elements  by  evaluating  partial 
derivatives  as  described  above.  The  sensitivity  of  the  optimal  solution 
can  also  be  compared  to  the  coefficients  of  the  empirically  derived  sys¬ 
tem  evaluation  models  in  order  to  identify  any  areas  where  additional 
data  may  be  required  to  improve  the  confidence  in  some  of  the  models. 

In  summary,  this  model  can  be  a  versatile  tool  for  a  variety  of 
system  analysis  techniques,  depending  upon  the  needs  of  the  user. 
Examples  of  several  of  these  will  be  presented  in  the  section  on 
applications. 


SECTION  III 


MODELING 


3.1  GENERAL 

In  this  section,  the  characteristics  of  the  dynamic  elements  will 
be  described  by  their  appropriate  differential  equations  or  transfer 
functions  in  terms  of  the  most  descriptive  and  useful  input  and  output 
parameters.  The  system  evaluation  models  used  in  determining  the  mea¬ 
sures  of  effectiveness  will  be  described  as  algebraic  combinations  of 
significant  statistics  of  the  dynamic  models  and  other  system  config¬ 
urations  and  baseline  input  data. 

Some  of  the  cost  models  will  be  developed  in  conjunction  with  the 
system  dynamics  models  in  order  to  preserve  continuity  in  the  present¬ 
ation.  These  will  be  referenced  in  the  development  of  total  system 
cost. 

3.2  DYNAMIC  ELEMENTS 

Of  the  eight  basic  system  elements  identified  in  Section  I,  seven 
of  these  elements  -  sea  state,  breakwater,  crane  platform,  container 
ship,  lighterage,  crane,  and  cargo  -  have  dynamic  characteristics  that 
influence  the  system  performance.  The  transfer  functions  of  the  break¬ 
water  and  the  floating  vessel  can  be  derived  independently  of  the  other 
system  dynamics  because  there  are  no  predominant  cross  coupling  effects 
The  crane  and  cargo,  however,  will  be  analyzed  as  a  single  entity 
because  the  dynamic  interaction  between  the  stiffness  of  the  crane  and 
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the  mass  of  the  cargo  contributes  strongly  to  the  frequency  character¬ 
istics  of  the  system. 

3.2.1  Sea  State 

As  discussed  earlier,  a  reasonable  and  useful  representation  of 

the  sea  state  is  its  power  density  spectrum.  A  number  of  researchers 

have  developed  analytical  expressions  for  sea  spectra  to  approximate 
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actual  wave  records.'  '  For  the  purpose  of  this  study,  any  spectrum 

that  approximates  the  type  of  seas  that  are  representative  of  the  COTS 

environment  would  be  adequate.  Many  researchers  accept  the  Pierson- 

Moscowitz  spectrum  for  fully  developed  sea  as  correct,  and  recent 

studies  by  the  Naval  Coastal  Systems  Center  have  determined  that  it  is 

(45) 

representative  of  the  majority  of  the  potential  COTS  environments. 

It  is  defined  as(46) 

S ( cl> )  -  (ag2/w5)  exp  (-8  u>4/u>4) 

for  values  of  from  zero  to  infinity,  where 
ct  =  8.10  x  10"3 
8  *  .74 

g  =  acceleration  of  gravity 
=  g/u 

u  =  wind  speed 

Any  set  of  g,  u,  and  u>  can  be  used  with  consistent  units.  The  accel¬ 
eration  of  gravity,  g,  is  taken  to  be  32.2  ft/sec2,  then  the  spectrum 

$(«)  =  (8.4/u>5)  exp  C795’530/u4a)4). 
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Since  the  significant  wave  height, 

H$  =  4.0 

the  wind  speed,  u,  can  be  expressed  in  terms  of  the  significant  wave 
height,  H$,  and  the  spectrum  can  be  written: 

S(u)  =  (8.4/u5)  exp  C33,6/Hs2w4), 

where  H  is  in  feet, 
s 

The  sea  state  power  density  spectrum  can  now  be  completely  defined 
by  the  significant  wave  height,  Hs-  This  provides  an  effective  descrip¬ 
tive  parameter  for  sea  state  with  a  strong  physical  connotation  for  the 
user.  The  Pierson-Moscowitz  sea  spectra  for  two  values  of  significant 
wave  height  are  shown  in  Figure  3.1. 

3.2.2  Breakwater 

As  discussed  earlier,  there  are  a  variety  of  designs  of  breakwaters. 

The  tethered  float  breakwater  (TFB),  shown  in  Figure  3.2,  was  chosen  as 

a  viable  candidate  for  this  study  because  there  have  been  a  number  of 
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theoretical  analyses  of  its  performance,  ’  '  wave  tank  model 

(481  (491 

tests,  '  engineering  design  and  deployment  studies,  and  a  full 

scale  prototype  test  is  currently  being  conducted  outside  San  Diego 
Bay. It  should  be  emphasized  that  the  COTS  model  is  not  restricted 
to  the  use  of  the  TFB  models.  Any  breakwater  design  whose  energy  trans¬ 
mission  ratio  can  be  calculated,  measured,  or  estimated  can  be  input  to 
the  model . 


rrvmY^f 

FT 

-4s 


FREEBOARD 


ORAFT 


‘MOORING  LINE 


FLOAT 

TETHER 


BALLAST 


SECTION  A-A 


Figure  3.2.  Tethered  Float  Breakwater 
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The  theoretical  energy  tranmission  ratios  (ETR)  as  a  function  of 
frequency  obtained  from  Scripps  Institute  of  Oceanography,  (Appendix  C) 
characteristically  appear  as  shown  below. 


The  tail  of  the  curve  would  indicate  that  high  frequency  energy  is 
transmitted  as  well  as  low;  however,  in  practice,  these  frequencies  are 
well  above  any  significant  energy  found  in  the  sea  and  well  above  the 
response  characteristics  of  the  floating  vessels.  For  this  reason,  in 
practical  applications,  the  TFB's  energy  transmission  ratio  (ETR)  can 
be  modeled  as  a  simple  low  pass  filter,  described  by  a  cutoff  frequency, 
u  ,  and  a  high  frequency  attenuation,  a. 
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Of  the  designs  reviewed  (Appendix  C),  nine  were  chosen  with  high 
frequency  attenuations  in  the  neighborhood  of  .05.  These  were  recom¬ 
mended  by  Scripps  as  being  the  most  promising  designs.  The  negative 
slopes  of  the  curves  beyond  the  cutoff  frequency  are  also  approximately 
equal.  Cutoff  frequency  ranges  from  1.0  to  2.0  rad/sec,  with  construc¬ 
tion  costs  ranging  from  $5250  to  $300  per  foot  respectively.  A  log-log 
plot  of  cost  versus  cutoff  frequency  is  shown  in  Figure  3.3,  indicating 
that  a  power  function  provides  a  reasonable  representation  of  this 
relationship. 

The  analytical  forms  of  ETR  and  breakwater  cost  (Appendix  C)  to  be 
used  in  the  model  are: 


!1 .0  ;  a)  <_  wi 

1.0  -  2.07  (u)  wj);  aj i  <  a)  <_  to 
.05 

where  “  Mc  ” 

and 

COST  («c)  *  5.4  [^-)'3 
c 

3.2.3  Floating  Elements 

A  vessel  oscillating  in  a  fluid  medium  under  the  influence  of 
random  waves  is  a  relatively  complex  phenomenon  that  is  difficult  to 
accurately  describe  and  predict  by  analytical  means.  The  vessel  is 
subjected  to  a  variety  of  forces.  There  are  forces  due  to  the  inertial 
effects  of  the  water  that  moves  in  harmony  with  the  oscillating  vessel. 


;  wc  <  w 


dollars 

ft 
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This  effect  is  referred  to  as  the  added  mass.  There  are  damping  forces 
that  remove  energy  from  the  vessel  through  secondary  wave  generation 
and  viscous  effects.  There  are  hydrostatic  restoring  forces  of  the 
water  pressure  acting  on  the  hull  of  the  ship,  and  there  are  forces 
resulting  from  the  action  of  the  incident  waves.  Consider  the  equation 
of  motion  for  vertical  heave: 


Mz  =  Fa  +  Fd  +  Fh  +  Fw 


where  z  «  vertical  displacement 
M  =  mass  of  the  ship 
F  *  added  mass  inertial  forces 

a 

Fd  =  damping  forces 

=  hydrostatic  restoring  forces 

F .  =  wave  forces . 
w 

In  general,  these  forces  are  not  expressible  as  linear  functions  of  z 
and  its  derivatives.  Observations,  however,  indicate  that  wall-sided 
vessels  in  moderate  sea  conditions  exhibit  linear  response  character¬ 
istics.  (24,25,26)  ,jynami-c  character  of  a  vessel  is  often  expressed 
in  terms  of  its  response  amplitude  operators  (RAO)  for  its  six  degrees 
of  freedom  -  heave,  surge,  sway,  putch,  roll,  and  yaw.  An  RAO  is  de¬ 
fined  as  the  ratio  of  the  amplitude  of  the  response  (either  translational 
or  angular)  to  the  amplitude  of  the  incident  waves  as  a  function  of 
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frequency.  A  ships's  RAO's  can  be  determined  experimentally  by  placing 
a  scale  model  in  a  wave  tank  and  recording  the  ship's  responses  to  known 
input  wave  conditions. 

Analytical  methods  have  been  developed  to  determine  ship  RAO's. 


The  most  completely  developed  and  frequently  used  technique  is  known  as 
strip  theory.  ’  ^  The  complex,  three-dimensional,  fluid  flow  problem 

of  a  six-degree-of- freedom,  oscillating  ship  is  approximated  by  a 
sumnation  of  several  simpler  two-dimensional  problems.  The  ship 
(Figure  3.4)  is  divided  into  a  number  of  transverse  slices,  each  rep¬ 
resented  by  a  two-dimensional  cylinder.  The  hydrostatic  and  dynamic 
forces  are  calculated  for  each  slice  independently  and  then  integrated 
along  the  length  of  the  vessel  to  obtain  the  total  effects  of  the 
forces  and  moments.  This  technique  has  been  Implemented  in  the  RELMO 
computer  program  developed  at  the  Navy's  Civil  Engineering  Lab¬ 
oratory.^®^ 

The  RELMO  program  was  originally  developed  to  predict  the  relative 
motions  between  points  on  two  ships.  In  order  to  do  this,  it  solves 
for  the  six  RAO's  simultaneously.  The  RAO’s  are  complex  functions  of 
frequency,  and  therefore,  contain  phase  as  well  as  amplitude  information. 
The  program  will  be  used  to  obtain  these  complex  RAO's  for  the  floating 
elements  in  the  analysis.  The  input  data  format  for  this  program  is 
given  in  Appendix  D. 

3.2.4  Crane/Cargo 


By  virtue  of  their  strong  Interaction,  the  crane  and  cargo  models 
cannot  be  developed  independently.  The  crane  is  subjected  to  a  variety 
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of  forces  during  the  offloading  process.  These  include  the  static 
weight  of  the  container,  acceleration  forces  -  some  of  which  are  amp¬ 
lified  by  the  structural  resonance  of  the  crane/cargo  system,  and  pen- 
dulation  loads  which  appear  primarily  as  side  loads  on  the  crane  boom. 
The  required  radius  of  operation  of  the  crane  will  depend  on  the  dim¬ 
ensions  of  the  ships  involved  and  the  system  configuration.  A  simpli¬ 
fied  dynamic  model  of  the  crane/cargo  system  is  shown  below. 


The  equations  of  motion  for  its  degrees  of  freedom  are  coupled. 
Linearization  and  simplification  of  the  equations  results  in  two  models 
that  describe  the  important  characteristics  of  the  system,  a  vertical 
stiffness  model  and  a  pendulation  model. 
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where 


/ 


M. 


9 

Z 


mass  of  cargo  and  hook  blocks 

stiffness  of  hoist  cable 

damping  of  hoist  cable 

effective  mass  of  moving  crane  structure 

equivalent  stiffness  of  crane  structure 

damping  of  crane  structure 

pendulation  length  of  load  line 

vertical  input  motions 

horizontal  input  motions 

pendulation  angle 

vertical  displacement  of  the  load 
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All  physical  systems  have  mechanisms  of  energy  dissipation,  such  as 

sliding  friction,  viscous  damping,  structural  damping,  and  the  like. 

Although  these  various  damping  mechanisms  are  represented  by  somewhat 

different  mathematical  descriptions,  for  very  lightly  damped  systems 

such  as  a  crane,  the  linear  approximation  of  equivalent  viscous  damping 

is  a  reasonable  representation  of  reality. ^ ^ 

Since  k,  >  k0  and  M„  <  for  the  crawler  cranes  under  consider- 
I  Z  sc 

( 36 ) 

ation,  the  vertical  stiffness  model  can  be  further  simplified,  ' 


where 


M  =  equivalent  suspended  mass 

k  =  total  vertical  stiffness 

C  =  equivalent  linear  damping  coefficient 


uv  =  vertical  displacement  of  the  boom  tip  if  the 
crane  were  a  rigid  body 
z  *  vertical  displacement  of  the  cargo. 

The  input  motions,  u(t),  induced  by  the  crane  platform  motions  are  re¬ 
lated  to  the  cargo  motion  z(t)  by  the  differential  equation. 


Mz  +  Cz  +  kz  *  Cuy  +  kuy 


i 


i 


I 


51 


The  transfer  function,  resulting  from  the  Laplace  transformation  is: 

z(s)  _  Cs  +  k _ 

uv(s)  Ms2  +  Cs  +  k 

The  dynamic  load  on  the  crane  resulting  from  the  suspended  mass  and  its 
acceleration  is  expressed  by, 

F  =  Mz  . 

When  this  expression  is  combined  with  the  transfer  function,  z(s)/uy(s), 
the  transfer  function  between  uy(s)  and  F(s)  is  found  to  be 

F(s) 
uv(s) 

When  this  force,  F,  is  added  to  the  weight  of  the  container,  W,  the 
result  is  a  force  approximating  the  vertical  component  actually  experi¬ 
enced  by  the  crane. 

A  similar  analogy  can  be  made  for  the  pendulation  model. 


s2M  (Cs  +  k) 
Ms2  +  Cs  +  k 
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leads  to 


Mx  +  Cx  +M£X  =  C*  +M£  , 

l  l  l  uh  l  uh 


e(s)  .  1 


uh(s) 


l  \  s2  +  2coj  s  +  to  2 

n_  n 
P  P 


where 


VgTT"  =  natural  frequency  of  pendulation. 


?  =  C  =  damping  ratio. 

2Vtslq/l 

The  side  load  at  the  point  of  suspension  can  be  approximated  by: 


Side  Load  *  Hg  sin  (e). 


For  conventional  boom  cranes,  manufacturers  recommend  that  this 

load  should  not  exceed  a  fixed  percentage  of  the  rated  vertical  load, 
(52) 

usually  2%.  For  cranes  equipped  with  a  rider  block  tagline  system, 

(53) 

this  figure  can  be  increased  to  as  much  as  1%.  ' 

The  vertical  stiffness  and  damping  coefficients  required  for  the 
above  models  can  be  easily  measured  on  the  cranes  of  interest.  These 
models  constitute  the  crane/cargo  structural  models  shown  in  the  system 
synthesis  flow  chart.  The  output  of  these  models  is  an  equivalent  load 
imposed  on  the  crane  by  virtue  of  its  operation  in  the  dynamic  environ¬ 
ment.  This  load  along  with  the  maximum  reach  requirements  of  the 
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system  configuration  define  the  crane  cost  by  the  formula 

Crane  Cost  (in  $1000.00)  =  W  R 

17356'  -  25.72  R) 

where  M  =  Load  in  lbs. 

R  =  Radius  of  operation  in  ft. 

This  relationship  is  developed  in  Appendix  A.l. 

The  motion  compensation  model  is  also  described  in  detail  in 
Appendix  A. 2.  It  appears  as  an  intermediate  system  between  the  input 
motion,  uv(t),  and  the  vertical  structural  model. 


The  transfer  function  for  this  motion  compensation  system,  devel¬ 
oped  in  Appendix  A. 2  is, 

(i  +  ts  )  (  i  +  ^-  +  -4 ) 

-  U)  U>  * 

ZQ(S)  =  _ nh  nh _  , 

uv(s)  (1  +  ts)(1  +  +  jV)+  KG(s) 

nh  nh 
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where 


t 


hydraulic  pump  swashplate  time  constant  »  .3  sec. 

j  / 

hydraulic  pump/motor  system  resonance  *  10  'sec. 


5  =  damping  ratio  =  -5 

K  *  system  feedback  gain 

G(s)  ■  series  compensation. 

The  system  feedback  gain,  K,  is  a  design  parameter  that  affects 
both  the  performance  and  the  cost  of  implementation  of  the  motion  com¬ 
pensation  system.  The  low  frequency  motion  attenuation  is  approximately 
equal  to  .  Increasing  X  reduces  the  low  frequency  motions  of  the 
load  and  increases  the  effective  bandwidth  of  the  system,  however,  this 
also  increases  the  horsepower  consumption  of  the  system  and  makes  sta¬ 
bilization  more  difficult.  The  difference  between  the  input  and  output 
velocities  is  the  hoisting  velocity  that  must  be  supplied  to  the  system. 
The  significant  value  of  this  hoist  line  velocity  can  be  calculated  by 
the  spectral  analysis  techniques  described  earlier.  The  transfer  func¬ 
tion  between  the  hoist  line  velocity  and  the  input  motion,  also  devel¬ 
oped  in  Appendix  A. 2,  is: 


V(s)  ,  _ sKG( s_) _ z _ 

u  (s)  (1  +  t  s)  (1  +Mi  + «3H2)  +  K  G(S) 

v  nh  h 

The  system  horsepower  can  be  calculated  by, 

hp,  - 

s  550 


HP 


s 


where 


significant  horsepower 


W  *  weight  of  cargo  in  lbs 

V$  =  significant  hoist  speed  in  ft‘/sec. 

A  motion  compensation  cost  model  of  the  form. 

Motion  Compensation  Cost  *  A  +  B  HP$  results, 

where  A  =  fixed  cost  of  major  motion  compensating  hardware 
components  and  installation  on  crane. 

B  *  variable  cost  associated  with  required  horsepower. 

These  cost  factors  are  addressed  in  further  detail  in  Appendix  A. 2. 

Because  effective  crane  operation  at  sea  requires  the  continuous 
positioning  of  a  pendulus  load,  a  simple,  idealized,  positioning  problem 
(shown  below)  is  examined  in  detail  in  Appendix  B.l. 


x(0)  =  x(0)  *  x ( 0 )  =  0 
u(0)  =  0 

u(tf)  =  x ( tf )  =  xf 
x(tf)  =  x ( tf )  =  0 

i(t!Uv,*ax 


Find  tf,  such  that  tf  is  minimum  over  all  u(t). 


It  is  shown  that  the  minimum  time  required  to  transfer  a  pendulum  from 
rest  to  a  nearby  position  at  rest  is  approximately  proportional  to  the 
square  root  of  the  pendulum  length.  This  relationship  will  be  used  in 
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the  development  of  analytical  expressions  for  system  productivity  in 
Appendix  B.2. 

3.3  MEASURES  OF  EFFECTIVENESS 

The  system  evaluation  model  uses  input  data  in  conjunction  with 
the  outputs  from  the  system  dynamics  model  to  generate  estimates  of 
system  cost  and  productivity.  The  following  sections  describe  how  the 
models  of  the  various  system  elements  are  combined  to  produce  the  values 
for  these  system  measures  of  effectiveness. 

3.3.1  Cost 

The  total  system  cost  is  the  sum  of  the  costs  of  the  various  sub¬ 
systems.  This  includes  the  breakwater  cost,  basic  crane  cost,  add¬ 
itional  crane  costs  associated  with  the  riderblock  or  motion  compensat¬ 
ing  hardware,  and  the  system  configuration  cost. 

The  breakwater  cost  model  has  been  described  along  with  the  break¬ 
water  dynamic  model  in  Section  3.2.2.  The  output  of  the  breakwater 
cost  model  is  cost  per  unit  length  in  dollars  per  foot.  This  must  be 
multiplied  by  the  length  of  breakwater  required  for  the  particular  con¬ 
figuration  under  consideration  to  obtain  the  total,  installed,  break¬ 
water  cost.  This  cost  is  based  on  domestic  material  costs  and  civilian 
labor  rates  in  effect  during  1977. 

The  basic  crane  and  motion  compensation  cost  models  have  been  given 
in  Section  3.2.4  along  with  the  crane  dynamic  models.  The  load  used  in 
the  equation  is  the  larger  of  either  the  vertical  dynamic  load  or  the 
effective  pendulation  load.  The  vertical  dynamic  load  is  the  static 
weight  of  the  cargo  plus  the  maximum  expected  vertical  inertial  load. 


The  effective  pendulation  load  is  the  maximum  expected  side  load  divided 
by  the  allowable  percent  side  load  for  the  particular  crane  under  con¬ 
sideration,  (.02  for  standard  cranes  and  .07  for  cranes  equipped  with  a 
rider  block). 

The  cost  of  the  Rider  Block  Tagline  System  (RBTS)  is  assumed  to 
be  $25,000  plus  one  tenth  of  the  basic  crane  cost.  This  figure  was  ob¬ 
tained  from  the  costs  of  a  prototype  RBTS  that  was  installed  on  a  small 
crane  and  a  full  scale  RBTS  that  is  currently  being  manufactured  for  a 
much  larger  crane. 

The  motion  compensation  cost  is  assumed  to  be  a  fixed  cost  of 
$250,000  plus  $400  per  horsepower  required  by  the  motion  compensation 
system^  (Appendix  C). 

The  system  configuration  cost  represents  all  other  fixed  costs 
that  are  associated  with  the  basic  system  configuration.  These  include 
the  costs  of  ships,  moorings,  and  other  hardware  as  well  as  the  indirect 
costs,  if  applicable,  of  transportability,  availability,  and  reliability 
that  are  further  descriptive  of  the  configuration.  This  information  is 
assumed  to  be  a  constant  for  a  particular  configuration. 

3.3.2  Productivity 

Productivity  can  be  defined  as  the  number  of  containers  offloaded 
in  a  given  period  of  time,  containers  per  hour  or  containers  per  day 
being  the  most  commonly  used.  It  is  a  somewhat  less  well  defined  par¬ 
ameter  than  system  cost.  There  has  been  very  little  work  done  in  dev¬ 
eloping  functional  relationships  between  other  system  parameters  and 
productivity  with  the  exception  of  one  report  that  used  linear 


regression  techniques  to  fit  the  times  required  to  perform  various 
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portions  of  the  offload  cycle  to  relative  vertical  displacement. 

The  offload  cycle  time  is  difficult  to  define  analytically  because 
it  is  a  function  of  so  many  variables,  many  of  which  are  not  readily 
quantifiable.  A  list  of  the  most  important  of  these  include: 

a.  Operator  Proficiency 

b.  System  Configuration 

c.  Sea  State 

d.  Vertical  Relative  Motion 

e.  Lateral  Relative  Motion 

f.  Maximum  Crane  Speeds 

g.  Pendulum  Length 

h.  Weight  of  the  Cargo 

i.  Load  Placement  Accuracy  Requirements 

j.  Communication  to  the  Crane  Operator 

k.  Weather  Conditions  other  than  Wave  Action 

Only  two  military  exercises  have  been  conducted  in  which  any 
significant  amount  of  productivity  data  was  collected.  These  were  the 
Over-the-Shore  Discharge  of  Cargo  II  (OSDOC  II)  exercises ^  conducted 
in  October  of  1972  and  the  Loglstics-Over-the-Shore  (LOTS)  exercises^16) 
conducted  in  August  of  1977.  During  OSDOC  II,  container  handling  oper¬ 
ations  were  performed  in  sea  states  1  and  2.  Data  were  taken  for  a 
variety  of  configurations.  It  was  observed  that  cranes  with  shorter 
booms  were  more  productive  than  cranes  with  longer  booms.  The  longer 
booms  resulted  in  longer  pendulums  in  most  cases.  In  both  exercises. 
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the  significant  time  losses  during  a  cycle  were  in  positioning  of  the 
load  over  the  ship  or  the  lighter.  These  observations  lead  to  the 
hypothesis  of  a  cycle  time  model  of  the  form: 

CT  "  Tm  +  Ts  +  T£  * 

where  CT  *  Total  Cycle  Time 

Tm  *  Minimum  practical  time  to  move 
a  container  in  complete  calm. 

Ts  »  Average  time  lost  over  the  ship. 

=  Average  time  lost  over  the  lighter. 

After  the  OSDOC  II  exercises,  several  observers,  knowledgable 
in  the  field  of  marine  operations  and  Naval  material  handling,  extrap¬ 
olated  the  data  to  generate  estimates  of  offload  cycle  times  for  sea 
(551  (561 

states  1  through  4.  '  The  minimum  time,  T  '  was  subtracted  from 

these  estimates  to  obtain  the  total  lost  time.  A  review  of  the  data 
indicated  that  the  time  lost  over  the  ship,  T$,  was  approximately  70% 
greater  than  the  time  lost  over  the  lighter,  T 

These  lost  times  were  hypothesized  to  have  a  form: 

T  -  A  aalB 

where 

A  *  constant 

o  a  relative  vertical  displacement  between 
the  cargo  and  the  ship  or  lighter 
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l  *  pendulum  length 
a&S  *  exponents 

Vertical  displacement  is  not  the  only  motion  that  affects  the  con¬ 
tainer  transfer  rate,  but  a  high  correlation  has  been  shown  to  exist 
(541 

between  the  two  and  the  other  motions  tend  to  increase  along  with 
vertical  displacement. 

The  pendulum  length  is  closely  related  to  the  control ability  of 
the  load.  Each  container  positioning  operation  is  actually  the  result 
of  a  finite  number  of  small  repositionings  and  minor  adjustments.  The 
number  required  increases  with  the  amount  of  relative  motion.  It  is 
assumed  that  an  experienced  crane  operator  might  perform  these  position¬ 
ing  tasks  in  some  suboptimal  manner  that  is  approximately  proportional 
to  the  optimal  solution  as  derived  in  Appendix  B.l.  Therefore,  the 
exponent,  B,  in  the  lost  time  expressions  is  chosen  to  be  1/2. 

A  typical  container  offloading  configuration  was  chosen  and  the 
relative  vertical  displacements  were  obtained  from  the  RELMO  computer 
program  for  the  various  sea  states.  This  information  is  plotted  against 
the  lost  time  data  in  Appendix  B.2  to  obtain  the  remaining  variables  in 
the  lost  time  expressions.  The  resulting  expressions  are: 


and 


Ts  •  .82asV 


JZ  = 


The  system  productivity  is  measured  in  containers  per  20  hour  day. 


Four  hours  are  allowed  each  day  for  maintenance  of  the  equipment. 
Therefore, 


0LK=^Oa 


where. 


OLR  =  offload  rate 
n  *  efficiency  factor. 


The  efficiency  factor  accounts  for  the  time  lost  moving  hatch- 
covers,  repositioning  the  cranes,  mooring  the  lighters  and  the  like. 
Preliminary  indications  from  the  LOTS  test  indicate  that  this  effi¬ 
ciency  factor  was  approximately  65%. 

3.4  COMPUTER  PROGRAM 

A  computer  program  has  been  written  (Appendix  E)  to  implement  the 
models  that  have  been  developed.  It  is  written  in  Fortran  and  makes 
extensive  use  of  the  subroutine  structure  to  permit  the  substitution  of 
alternative  models  for  the  various  system  elements  and  to  facilitate 
the  reconstruction  of  the  output  for  a  variety  of  applications. 

The  ship  response  operators  are  calculated  by  a  separate  program 
and  stored  on  tape  or  disk  files.  This  data  is  called  and  read  by  the 
main  program  at  the  beginning  of  each  run.  This  Is  done  because  the 
evaluation  of  the  response  amplitude  operators  for  the  six  degrees  of 
freedom  of  a  floating  vessel  Is  significantly  more  expensive  than  run¬ 
ning  the  remainder  of  the  system  dynamics  and  system  evaluation  models. 
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The  subroutine  names  and  their  functions  in  the  main  program  are: 

RAOG  -  reads  RAO's  for  the  required  number  of  vessels. 

PIERMOS  -  Calculates  PIERSON-MOSCOWITZ  Sea  Spectrum 
for  given  significant  wave  height 
TFB  -  calculates  Energy  Transmission  Ratio  and  cost 
for  Tethered  Float  Breakwater 
CRANE  -  calculates  structural  and  pendulation  transfer 
functions  for  crane 

MOCOMP  -  calculates  motion  compensation  transfer 
function 

CRCOST  -  calculates  total  crane  cost 
PROD  -  calculates  system  productivity  in  offload  rate 
INTEG  -  calculates  area  under  a  power  density  spectrum 
EVAL  -  combines  the  input  data  with  the  other  subroutines  to 
generate  total  crane  cost,  breakwater  cost,  total 
system  cost,  and  offload  rate. 

Input  data  can  be  provided  through  two  namelists  -  CONFIG  and 
VARIAB.  The  system  configuration  data  contained  in  CONFIG  includes: 

XCG  -  Relative  position  of  the  ship  and  lighter  with  respect 
YCG  to  the  crane  platform 
XR  -  Points  of  Interest  on  the  ship  and  the 
YR  lighter 
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CFCOST  -  Fixed  configuration  cost 

BWL  -  Breakwater  length 

MCO  -  Motion  compensation  option 

0  »  standard  crane 

1  *  rider  block  tagline  system 

2  *  rider  block  with  motion  compensation 
NC  -  Number  of  Cranes 

PMT  -  Practical  minimum  cycle  time 
EFF  -  Efficiency  factor 

The  other  variables  that  are  contained  in  VARIAB  are: 

WC  -  Container  Weight 

VS  -  Vertical  Stiffness  of  Crane  Structures 
FK  -  Motion  compensation  feedback  gain 
BCF  -  Breakwater  cutoff  frequency 
H13  -  Significant  wave  height 

The  program  has  been  run  on  a  commercial  time  sharing  system 
(Control  Data  Corporation's  Cyber  174).  The  program  can  be  compiled 
and  180  variations  of  a  basic  system  configuration  can  be  evaluated  for 
approximately  $24.00.  The  RAO's  for  the  six  degrees  of  freedom  of  a 
floating  vessel  can  be  calculated  and  stored  on  a  file  for  approximately 
$30.00. 


SECTION  IV 
APPLICATIONS 

4.1  GENERAL. 

The  COTS  computer  model  developed  in  the  preceding  sections  can  be 
used  in  a  variety  of  ways.  In  this  section,  a  number  of  potential  appli¬ 
cations  will  be  demonstrated  using  the  model.  This  is  followed  by  two 
actual  applications  using  test  data.  A  comparative  analysis  between  the 
model  and  the  actual  test  data  then  is  presented. 

4.2  THEORETICAL  DEMONSTRATIONS 

For  the  purpose  of  demonstrating  potential  applications  of  the  model, 
an  example  system  configuration  will  be  investigated.  It  is  shown  in 
Figure  4.1.  "he  crane  platform  is  a  C7  class  container  ship,  heavily  loaded 
for  a  lower  roll  resonance  providing  better  stability  in  the  dynamic  ocean 
environment.  The  container  ship  is  a  C4  class,  moderately  loaded,  simulat¬ 
ing  a  partially  unloaded  vessel.  The  lighter  is  an  LCM-8,  a  fairly  light 
landing  craft  capable  of  carrying  one  container  at  a  time. 

The  ships  are  in  quartering  seas.  That  is,  the  direction  of  the 
waves  is  135°  counter  clockwise,  measured  from  the  stern  of  the  vessels. 

The  center  of  gravity  of  the  crane  platform  is  used  as  a  reference 
for  the  definition  of  the  system  configuration.  The  center  of  gravity 
of  the  container  ship  and  the  lighter  is  located  with  respect  to  the 
right  handed  coordinate  system  whose  origin  is  the  center  of  gravity  of 
the  crane  platform.  The  center  of  rotation  of  the  crane  is  also  located 
with  respect  to  this  coordinate  system.  The  X  axis  faces  forward  along 
the  longitudinal  axis  of  the  ship.  The  X-Y  plane  is  parallel  to  the 
horizon,  with  the  Z  axis  pointing  vertically  upward.  The  X  and  Y 
coordinates  of  the  container  ship  are  0  and  -110  feet  (0,-33. 5m) 
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Figure  4.1.  An  Example  COTS  Configuration 
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respectively.  The  coordinates  of  the  lighter  are  60  and  70  feet  (18.3, 
21.3  m).  The  center  of  rotation  of  the  crane  is  at  80  and  0  feet  (24.4, 
0  m).  The  reference  point  of  the  containership  and  the  lighter  used  in 
the  estimation  of  reach  and  motion  requirements  is  defined  from  sec¬ 
ondary  coordinate  systems  whose  origins  are  at  the  center  of  gravity  of 
the  respective  vessels.  The  X  and  Y  coordinates  of  the  points  of  the 
container  ship  and  the  lighter  are  120  and  -20  feet  (36.6,  6.1  m)  and 
0  and  0  feet  (0,  0  m)  respectively. 

The  other  input  variables  used  to  describe  this  configuration  are: 

Fixed  Cost  (CFCOST)  =  $1,000,000 
Breakwater  Length  (BWL)  *  1200  feet  (365.8  m) 

Motion  Compensation  Option  (MCO)  ■  0  *  (standard  crane) 

Number  of  Cranes  (NC)  =  1 

Practical  Minimum  Time  (PMT)  =4.0  minutes 

Efficiency  (EFF)  =  0.65 

Container  Weight  (WC)  =  40,000  lbs  (18,144  kg) 

Vertical  Stiffness  (VS)  =  3500  Ibs/in  (395  N/m) 

Motion  Compensation  Gain  (FK)  =  0 

Breakwater  Cutoff  Frequency  (BCF)  =  4.0  rad/sec 

Significant  Wave  Height  ( HI 3)  =  6.0  feet  (1.83  m) 

These  data  along  with  the  configuration  coordinates  and  the  re¬ 
sponse  amplitude  operators  for  the  three  vessels  constitute  the  initial 
input  data  for  the  program.  These  variables  can  be  altered  within  the 
program  in  order  to  evaluate  a  variety  of  alternative  system  variations. 
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4.2.1  EFFECT  OF  SEA  STATE 

A  useful  application  of  the  model  is  to  predict  container  offload 
rates  of  a  given  system  configuration  for  a  range  of  sea  states.  This 
can  easily  be  done  by  incrementing  the  significant  wave  height  and 
solving  for  the  offload  rate.  It  is  also  possible  to  repeat  this  ex¬ 
ercise  for  other  configurations. 

Motion  compensation  options  (MCO)  of  0  and  1,  which  is  a  standard 
crane  and  a  crane  with  a  rider  block  tagline  respectively,  are  used, 
along  with  significant  wave  heights  ( HI 3 )  of  1  to  8  feet.  The  initial 
condition  of  BCF  =  4.0  rad/sec  is  interpreted  by  the  TFB  subroutine  as 
no  breakwater,  since  this  cutoff  frequency  is  sufficiently  high  to  have 
no  effect. 

The  output  of  this  computer  run  is  given  in  Table  4.1.  Figure  4.2 
is  a  plot  of  these  results.  Container  offload  rate  is  shown  versus  wave 
height  for  a  standard  crane  and  a  crane  with  a  rider  block  (RBTS).  The 
five  to  six  foot  wave  heights  represent  an  upper  sea  state  3,  the 
highest  operating  range  of  a  container  offload  and  transfer  system 
using  the  LCM-8  lighters.  The  seaworthiness  of  these  lighters  becomes 
an  operational,  limiting  factor,  although  the  theoretical  offload  rates 
are  shown  above  sea  state  3  for  comparison  purposes.  Container  offload 
rate  decreases  rapidly  with  increasing  wave  height  for  wave  heights 
above  about  one  foot  for  both  standard  and  RBTS  cranes.  In  a  sea  state 
3,  the  offload  rate  of  a  standard  crane  is  reduced  to  about  25%  of  its 
value  in  calm  water.  As  sea  conditions  worsen,  the  offload  rate  of  the 
standard  crane  decreases  more  rapidly  than  the  RBTS  crane. 
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2.  Effect  of  Wave  Height  on  Productivity  of  Exampl 
COTS  with  Standard  and  RESTS  Equipped  Cranes 
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Figure  4.3  is  a  plot  of  the  percentage  increase  in  productivity  and 
cost  of  the  RBTS  crane  over  the  standard  crane  versus  wave  height. 

The  benefits  of  the  RBTS  crane,  in  terms  of  increased  productivity, 
increases  with  increasing  wave  height.  In  a  sea  state  3,  the  RBTS 
crane  has  an  offload  rate  that  is  approximately  50%  greater  than  a 
standard  crane. 

The  added  cost  of  the  RBTS  equipment  is  approximately  12%  of  the 
standard  crane  cost.  This  results  in  a  2-1/2  foot  wave  height  as  the 
economic  break  even  point  in  the  decision  to  make  use  of  an  RBTS  crane. 
At  wave  heights  less  than  this,  the  increased  cost  of  the  RBTS  hardware 
over  the  standard  crane  is  greater  than  the  percentage  increase  in 
productivity  derived  from  its  use. 

4.2.2  EFFECT  OF  BREAKWATERS 

Another  application  of  the  COTS  program  is  to  examine  the  effect 
of  variations  in  breakwater  design  on  system  productivity.  This  is 
easily  done  by  varying  the  breakwater  cutoff  frequency  (BCF)  for  a 
given  sea  condition.  The  changes  to  the  computer  program  required  to 
implement  this  are  given  in  Appendix  E. 

Once  again,  the  standard  crane  and  the  rider  block  option  are 
compared.  The  breakwater  cutoff  frequency  (BCF)  is  varied  from  4.0 
rad/sec  to  0.5  rad/sec  in  increments  of  0.5  rad/sec.  The  incident 
significant  wave  height  (H13)  in  all  cases  is  the  initial  input  value 
of  6.0  feet.  The  output  of  this  computer  run  is  given  in  Table  4.2. 
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Figure  4.3.  Percent  Increase  in  System  Cost  and  Productivity 
of  RBTS  Over  Standard  Cranes  Versus  Wave  Height 
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Figure  4.4  Is  a  plot  of  offload  rate  for  the  standard  and  RBTS 
cranes  versus  breakwater  cutoff  frequency  from  these  results.  Once 
again,  as  in  the  previous  example,  the  RBTS  crane  is  seen  to  provide  a 
larger  offload  rate  than  the  standard  crane  for  the  same  conditions. 

As  the  breakwater  cutoff  frequency  is  reduced,  thereby  reducing  the 
amount  of  wave  activity,  the  offload  rate  of  both  crane  designs  in¬ 
creases.  This  increase  is  very  slight  between  4  rad/sec  and  1.5  rad/sec 
but  rises  very  sharply  between  1.5  rad/sec  and  .5  rad/sec.  This  result 
is  not  too  surprising  upon  review  of  the  sea  spectrum  for  six  foot  waves 
(Figure  3.1).  The  majority  of  its  energy  is  found  at  frequencies  below 
1.5  rad/sec.  At  1.5  rad/sec,  the  productivity  of  both  cranes  has  in¬ 
creases  less  than  8 %  over  the  value  for  no  breakwater  (4.0  rad/sec). 

At  1.0  rad/sec,  however,  the  standard  crane  shows  a  77%  increase  and 
the  RBTS  crane,  a  62%  increase. 

These  productivity  figures  are  accompanied  by  costs  that  also  in¬ 
crease  most  rapidly  where  the  greatest  benefits  are  derived.  (See 
Table  4.2).  At  first  glance,  these  costs  appear  to  make  the  breakwater 
far  too  expensive  in  comparison  to  its  effect  on  the  system.  It  should 
be  noted  that  a  single  breakwater  can  shelter  more  than  one  crane.  Two 
cranes  could  easily  operate  from  the  same  crane  platform,  offloading 
the  same  containershlp.  As  the  breakwater  cost  increases,  the  crane 
cost  is  reduced,  because  the  reduction  in  the  transmitted  wave  height 
reduces  the  dynamic  loads  experienced  by  the  crane.  This  savings,  times 
the  number  of  cranes,  can  be  deducted  from  the  cost  of  the  breakwater  in 
making  comparisons. 
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Figure  4.4.  Effect  of  Breakwater  Cutoff  Frequency  on  Productivity 
of  Example  COTS  with  Standard  and  RBTS  Cranes 
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The  use  of  a  1.0  rad/sec  breakwater  in  the  example  with  two 
standard  cranes  results  in  an  increase  of  53%  in  the  total  system  cost. 
This  is  more  than  offset  by  the  increased  productivity  of  77%.  On  the 
other  hand,  the  selection  of  a  2.0  rad/sec  breakwater  for  this  system 
would  have  resulted  in  a  total  system  cost  increase  of  over  3%  with  a 
corresponding  increase  in  productivity  of  less  than  1%. 

Therefore,  it  appears,  that  if  breakwaters  are  to  be  used  in  such 
a  system,  that  one  large  enough  to  have  a  sufficiently  low  cutoff  fre¬ 
quency  would  be  more  cost  effective  than  the  smaller,  less  costlier 
designs . 

4.2.3  TRADEOFF  ANALYSIS 

Another  use  of  the  COTS  computer  model  is  in  evaluating  tradeoffs 
between  two  or  more  competing  technologies.  For  the  example  COTS  con¬ 
figuration  under  investigation,  the  tradeoff  between  motion  compensating 
cranes  and  breakwaters  will  be  examined. 

For  a  two  parameter  system  such  as  this,  the  brute  force  approach 
of  evaluating  the  system  over  the  ranges  of  both  independent  variables 
provides  the  designer  with  a  quick  look  at  the  nature  of  the  tradeoffs 
and  the  approximate  values  and  locations  of  any  optimal  solutions. 

The  required  lines  of  computer  code  to  implement  this  two  para¬ 
meter  evaluation  technique  are  listed  in  Appendix  E.  The  breakwater 
cutoff  frequency  is  varied  from  0.8  to  4.0  rad/sec.  The  motion  com¬ 
pensation  gain  is  varied  from  0.0  to  10.  The  resulting  offload  rate, 
total  system  cost,  and  cost  effectiveness  (total  system  cost/offload 
rate)  are  shown  in  Figures  4.5,  4.6,  and  4.7  respectively. 
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SYSTEM  COST/ PRODUCTIVITY 
( $1 OOO/BOX/DAY ) 

130  j 


Figure  joined  Effect  of  Motion  Compensation  Gain  and  Breakwater 
Cutoff  Frequency  on  Cost  Effectiveness  of  Example  COTS 
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The  effect  of  breakwater  cutoff  frequency  on  container  offload  rate 
that  was  observed  in  the  previous  example  (Section  4.3)  is  seen,  in 
Figure  4.5,  to  be  characteristic  of  the  effect  of  breakwaters  on  offload 
rate  regardless  of  the  presence  or  magnitude  of  motion  compensation 
gain.  This  characteristic  is  also  observed  with  respect  to  system  cost 
as  shown  in  Figure  4.6. 

The  effect  of  motion  compensation  gain  on  offload  rate  is  quite 
different  from  that  of  breakwaters.  The  addition  of  a  small  amount  of 
gain  results  in  a  pronounced  increase  in  offload  rate.  Further  in¬ 
creases  in  gain  result  in  successively  smaller  increases  in  offload 
rate.  This  is  explained  by  the  manner  in  which  gain  affects  the  motion 
compensating  system's  ability  to  track  the  motions  of  the  crane  plat¬ 
form.  A  measure  of  the  control  system's  capability  is  its  steady  state 
error  to  a  step  input  which  equals  1/K+l ,  where  K  is  gain.  An  increase 
in  gain  from  0  to  3  results  in  a  reduction  in  this  error  from  100%  to 
25%  or  75%.  An  additional  increase  in  gain  to  9,  only  results  in  an 
additional  15%  reduction  in  steady  state  error. 

Figure  4.6  shows  that  the  effect  of  motion  compensation  gain  on 
system  cost  is  not  significant.  Once  the  initial  cost  of  the  motion 
compensating  hardware  is  absorbed  in  the  total  system  cost,  the  in¬ 
creased  cost  associated  with  the  additional  horsepower  requirements, 
necessitated  by  increased  gain,  is  not  a  major  factor. 

Figure  4.7  provides  an  interesting  overview  of  the  combined  effects 
of  breakwater  cutoff  frequency  and  motion  compensation  gain  on  the  cost 
effectiveness  of  the  system.  The  cost  per  box  per  day  increases  with 


decreasing  breakwater  cutoff  frequency  between  4.0  rad/sec  and  1.5  rad/ 
sec.  At  this  point  the  cost  effectiveness,  as  measured  in  dollars  per 
box  per  day,  falls  sharply  with  decreasing  cutoff  frequency  to  approx¬ 
imately  1.0  rad/sec.  Decreasing  cutoff  frequencies  below  1.0  rad/sec 
result  in  rapidly  increasing  costs  per  unit  offload  rate. 

The  effect  of  motion  compensation  gain  on  cost  effectiveness  is 
monotonic.  At  all  points,  an  increase  in  gain  results  in  a  reduction 
in  cost  per  unit  offload  rate,  although  the  predominant  reductions 
occur  at  the  lower  values  of  gain. 

These  results  would  lead  a  COTS  designer  to  utilize  the  maximum 
motion  compensation  capability  at  his  disposal  and  fill  in  with  break¬ 
water  as  required  to  obtain  the  necessary  offloading  capability.  Two 
basic  COTS  configurations  appear  as  potentially  viable  candidates  from 
minimization  of  the  system  cost  effectiveness: 

a)  Maximum  motion  compensation  gain  and  no  breakwater,  and 

b)  Maximum  motion  compensation  gain  and  a  breakwater  with 
a  cutoff  frequency  in  the  neighborhood  of  1.0  rad/sec. 

System  (a)  provides  an  offload  rate  of  132  boxes/day  at  a  cost  of 
13.8  million  dollars  for  a  six  month  duration.  System  (b)  provides  an 
offload  rate  of  208  boxes/day  at  a  cost  of  20.7  million  dollars. 

The  two  parameter  system,  as  shown  above,  can  be  evaluated,  the 
resulting  surface  plotted,  and  system  tradeoffs  readily  identified. 

This  approach  breaks  down,  however,  for  systems  with  three  or  more 
variables.  In  these  cases,  optimal  solutions  could  be  found  by  gradient 
search  techniques  such  as  a  steepest  descent  algorithm. 
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addition  of  penalty  functions  to  the  measure  of  effectiveness,  permits 
the  use  of  this  technique  for  constrained  optimization. 

A  simple  minimization  routine  is  performed  on  the  cost  effective¬ 
ness  function.  Partial  derivatives  of  the  cost  effectiveness  with  re¬ 
spect  to  breakwater  cutoff  frequency  and  motion  compensation  gain  are 
numerically  evaluated  by  small  variations  about  their  present  values. 
These  numerical  partial  derivatives  are  used  to  calculate  the  unit 
vector  in  the  direction  of  the  gradient.  A  fixed  size  step  is  taken  in 
the  direction  of  the  gradient.  If  a  reduction  in  cost  effectiveness 
is  found,  the  process  is  repeated.  If  a  reduction  is  not  observed,  the 
step  size  is  repeatedly  cut  in  half  until  a  reduction  is  observed.  If 
the  step  size  is  reduced  to  a  sufficiently  small  value  with  no  observed 
reduction  in  cost  effectiveness,  the  process  is  stopped  and  the  coor¬ 
dinates  and  related  system  parameters  are  printed.  The  required  lines 
of  computer  code  are  listed  in  Appendix  E.  Depending  upon  the  initial 
starting  points,  the  routine  either  converges  to  the  local  minimum  at 
a  cutoff  frequency  of  4.0  rad/sec  and  a  gain  of  10.0  or  the  global  min¬ 
imum  where  cutoff  frequency  equals  1.05  rad/sec  and  gain  equals  10.0. 

An  example  of  the  use  of  the  COTS  model  in  a  constrained  minimiza¬ 
tion  problem  is  to  find  the  least  cost  configuration  that  will  provide 
an  offload  rate  of  at  least  170  boxes/day.  This  is  done  by  adding  to 
the  cost  function,  a  penalty  function  of  the  form: 
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where  OLR  is  the  offload  rate.  This  penalty  function  becomes  very  large 
for  offload  rates  less  than  170  boxes/day,  but  almost  vanishes  for  off¬ 
load  rates  greater  than  170  boxes/day. 

The  lines  of  computer  code  required  to  execute  this  procedure  are 
also  given  in  Appendix  E.  The  solution  obtained  for  this  example  is: 

Motion  Compensation  Gain  =  10.0 
Breakwater  Cutoff  Frequency  =  1.196  rad/sec 
Offload  Rate  =  170.9  boxes/day 
System  Cost  =  $17.3  million 

4.3  COMPARATIVE  ANALYSIS  WITH  ACTUAL  TEST  DATA. 

Program  COTS  was  tested  for  validity  by  comparative  analysis  with 
actual  data  from  two  COTS  test  series:  the  RTBS  harbor  tests^^ ,  and  the 
TCDF  ship  tests  (TECHEVAL  and  0PEVAL).(60’  61  *  62) 

4.3.1  HARBOR  TESTS  OF  CRANE  WITH  RBTS. 

With  the  Manitowoc  4100W  COTS  prototype  crane  installed  on  an  Army 
B-DeLong  Barge,  container  handling  tests  were  conducted  in  Port  Hueneme,  CA, 
in  November  and  December  1979.  Tests  were  first  conducted  at  dockside, 
where  there  was  no  appreciable  platform  motion.  Data  from  the  dockside 
evolutions  provided  an  estimated  PMT  of  3.7  minutes  for  input  to  the  COTS 
computer  program. 

Barge  RAOs  were  computed  with  program  "RAOS".  Some  estimates  were 
required  for  section  geometries  needed  by  the  program.  The  efficiency  factor 
for  transfer  was  assumed  to  be  1.0  since  the  test  was  fully  controlled. 

Following  the  dockside  tests,  the  barge  was  moved  to  the  inner  harbor. 
Figure  4.8  shows  the  transfer  scenario.  Five  transfer  sequences  were  conducted: 
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three  with  RBTS  and  two  without.  Table  4.3  contains  important  configuration 
parameters  used  in  the  computer  programs. 

A  review  of  the  test  results  revealed  an  anomaly  which  skewed  the  overall 
transfer  rate  average  for  the  RBTS  operation.  Two  of  the  RBTS  sequences  were 
performed  by  crane  operator  "A"  and  the  other  by  "B",  while  each  performed  one 
of  the  without-RBTS  sequences.  Although  B's  non-RBTS  sequence  was  significantly 
"better"  than  A's  (13.73  containers/hour  compared  with  11.67  for  A),  his  RBTS 
sequence  was  far  less  productive  than  A's  (9.24  compared  with  11.24  and  13.22 
for  A).  With  sufficient  operator  training,  even  at  low  seastates,  productivity 
can  be  expected  to  increase  when  RBTS  is  used.  Therefore,  it  can  only  be  assumed 
that  Operator  B  was  not  ready  for  a  record- purposes  sequence  using  the  RBTS. 
Meanwhile,  Operator  A's  RBTS  performance  of  11.24  and  13.22  (average  12.15)  was 
better  overall  than  his  sequence  without  RBTS  (11.67).  This  was  to  be  expected, 
given  the  designed  contribution  of  RBTS  to  the  system,  and  providing  that  the 
operator  is  sufficiently  trained  and  ready  for  the  RBTS  operation. 

It  is  recognized  that  the  sample-size  is  too  small  for  any  conclusions  to 
be  drawn  on  the  validity  of  the  computer  program  in  this  test.  In  view  of  the 
above,  however,  it  was  considered  appropriate  to  reduce  the  sample  even  further 
by  using  only  Operator  A's  performance  for  the  analysis. 

Table  4.4  shows  a  comparison  between  Operator  A's  test  rates  and  the 
computer  simulation  rates.  Since  tests  were  conducted  only  in  wave  heights  of 
<1  foot,  a  test  rate  slope  cannot  be  obtained  for  comparison  with  the  simulation 
rates  slope. 

Referring  back  to  Table  4.3  and  Figure  4.8,  the  "containership"  and  TCDF 
platform  are  barges.  Together,  they  produce  significantly  greater  movement,  as 
wave  heights  increase,  than  do  the  ship  platforms  used  for  the  OPEVAL  (4.3.2., 
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TABLE  4.3.  RBTS  TEST  PARAMETERS 


CRANE  SHIP 

CONTAINER  SHIP 

LIGHTER 

Length,  ft 

150 

72 

90 

Beam,  ft 

60 

25 

30 

Draft,  ft 

3.5 

1.5 

2 

Displacement,  long  tons 

900 

77 

154 

Roll  period 

6.2 

2.5 

2.5 

CRANE 

Boom  length,  ft  200 
Vertical  stiffness  (VS),  lb/in.  3500 
Practical  minimum  transfer  (PMT),  sec  3.7 


TABLE  4.4.  RBTS  TEST  TRANSFER  RATE 
CONTAINERS  PER  HOUR 


WITHOUT  RBTS 

TEST 

SIMULATION 

WAVE  HEIGHT 

RATE 

RATE 

H  1/3,  Ft 

1 

11.67 

11.22 

2 

- 

2.95 

3 

- 

.94 

WAVE  HEIGHT 

WITH  RBTS 

TEST 

RATE 

SIMULATION 

RATE 

H  1/3,  Ft 

1 

12.15 

15.31 

2 

- 

4.82 

3 

- 

1.70 
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below).  Therefore,  the  Table  4.4  simulation  rates  change  dramatically  as  wave 
heights  increase.  Further,  Wave  Height  1  effects  the  ships  relatively  little, 
while  barge  movement  is  significantly  greater  than  the  ships  at  that  wave  height. 
This  explains  the  greater  improvement  using  RBTS  in  the  barge  scenario.  Other 
features  of  Table  4.4  worthy  of  note:  (1)  Wave  Height  1  rates  are  much  higher  than 
OPEVAL  rates  (Table  4.6),  primarily  due  to  the  use  of  different  efficiency  fac¬ 
tors;  (2)  the  "Without  RBTS"  test  rate  compares  favorably  with  the  simulation 
rate,  while  the  "With  RBTS"  test  rate  does  not,  being  almost  3  containers/hour 
lower.  Although  (as  stated  earlier)  the  small  sample  size  mitigates  against  con¬ 
clusive  findings,  this  could  indicate  high  potential  for  improvement  through 
RBTS  operator  training. 

4.3.2  TCDF  SHIP  TESTS  (TECHEVAL/OPEVAL) . 

These  tests  were  performed  in  the  Silver  Strand  Amphibious  Beaching  Area 


November-December,  1980.  The  crane  was  mounted  on  the  ex-LSD16  (ex-CABILDO) , 
the  containership  was  the  SS  LINCOLN  (a  C-4  ship)  and  the  lighter  was  modeled 
as  an  LCM-8.  Table  4.5  contains  the  configuration  parameters  and  Figure  4.9 


TABLE  4.5.  TCDF  TEST 

PARAMETERS 

CRANE  SHIP 

CONTAINER  SHIP 

LIGHTER 

Length,  ft 

458 

563 

51 

Beam,  ft 

72 

76 

21 

Draft,  ft 

17 

15 

2.5 

Displacement,  long  tons 

9000 

14000 

51 

Roll  Period 

8 

12.4 

3 

CRANE 

Boom  Length,  ft  200 
Vertical  stiffness,  lb/in.  3500 
Practical  minimum  transfer  time,  min.  5.2 
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shows  the  scenario.  The  PMT  was  estimated  by  considering  several  factors: 

(1)  time-motion  analysis  of  the  harbor  test  (4.3.1,  above);  (2)  the  physical 
lay-out  of  the  TCDF  scenario;  and  (3)  the  increased  demands  on  operator 
performance  imposed  by  the  open-sea  environment. 

The  tests  were  performed  in  sea  states  1  and  2  (wave  heights  H  1/3  13  ft). 
Efficiency  was  estimated  to  be  0.65  since  this  test  represented  actual  opera¬ 
ting  conditions.  Results  were  averaged  from  all  available  data,  presented  in 

Table  4.6,  and  graphed  on  Figure  4.10.  Simulation  data  from  the  "COTS"  computer 
program  also  is  provided  for  comparison. 

Figure  4.10  is  a  plot  of  "with  RBTS"  data  only,  since  the  tests  were 
conducted  using  the  RBTS  throughout.  The  slope  of  the  simulation  curve  compares 
favorably  with  overall  TECH/OPEVAL  results  (the  TECHEVAL /OP EVAL  curve). 
Comparative  slope  is  the  singlemost  important  factor  in  validation  of  the 
computer  program,  since  vertical  displacement  of  the  simulation  curve  is  almost 
solely  dependent  upon  selection  of  a  realistic  PMT  and  the  appropriate 
efficiency  factor. 

It  also  should  be  noted  that  in  low  sea  states  (such  as  in  these  tests) 
the  PMT  is  the  dominant  operational  factor  for  transfer  rate.  As  sea  state 
increases,  ship  motions  and  crane  parameters  become  more  contributive. 

It  is  recognized  that  sample  size  from  these  tests  is  too  small  for 
conclusive  statistical  analysis,  but  it  is  interesting  to  note  that  as  sample 
size  increases,  the  test  data  curve  approaches  the  slope  of  the  simulation 
curve.  This  is  demonstrated  in  Figure  4.10  by  plotting  two  test  data  curves: 
one  for  TECHEVAL  only  (128  samples),  and  one  for  TECHEVAL  and  OPEVAL  (324 
samples) . 

Comparative  analysis  between  the  test  results  and  the  simulation  rates. 
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when  all  of  the  above  factors  are  considered,  tends  to  validate  the  "COTS" 
computer  program.  However,  only  through  considerably  more  at-sea  experience 
and  the  availability  of  more  test  data  (making  possible  accurate  estimates  of 
PMT  and  efficiency  factor)  can  accurate  predictions  of  transfer  rates  be  made. 


SECTION  V 


CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  CONCLUSIONS 

A  methodology  has  been  developed  for  the  evaluation  of  the  cost 
and  effectiveness  of  a  complex  system  of  dynamic  elements  driven  by  a 
single  random  input.  It  has  been  shown  that  the  interactions  among  the 
system  elements,  their  environment,  and  the  system  measures  of  effect¬ 
iveness  can  be  grouped  into  two  types  of  relationships  -  system  dynam¬ 
ics  and  system  evaluation.  When  the  system  dynamics  are  describable  by 
linear  differential  equations  and  the  input  is  a  zero  mean,  Gaussian, 
random  process,  the  system  responses  are  completely  described  by  their 
power  density  spectra.  The  significant  and  maximum  expected  values  of 
those  responses  that  are  inputs  to  the  system  evaluation  relationships 
can  then  be  obtained.  These  evaluation  relationships,  that  define  the 
system  measures  of  effectiveness,  can  be  any  function,  linear  or  non¬ 
linear,  of  the  descriptive  parameters  of  the  system  configuration  and 
the  significant  and  maximum  expected  values  of  the  system's  dynamic  re¬ 
sponses  . 

In  particular,  this  methodology  was  developed  for  the  evaluation 
of  amphibious  material  handling  systems  and  applied  to  the  Navy's 
Container  Offloading  and  Transfer  system  (COTS).  An  example  system  con¬ 
figuration  has  been  examined  in  detail  in  order  to  demonstrate  the 
application  of  the  methodology  and  the  type  of  data  needed  to  develop 
the  necessary  system  relationships.  A  variety  of  applications  were 
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demonstrated,  including  determination  of  the  effects  of  sea  state  and 
breakwaters  and  analysis  of  the  tradeoff  between  breakwaters  and  motion 
compensating  cranes.  It  was  demonstrated  that  the  methodology  is 
easily  adaptable  to  multidimensional  optimization  by  the  use  of  numer¬ 
ical  techniques  such  as  the  steepest  descent  algorithm.  Further,  by  com 
oarative  analysis  with  a  limited  sample  of  actual  test  data,  it  was  dem 
onstrated  that  Program  "COTS"  is  a  potentially  useful  tool  for  predic¬ 
ting  crane  performance  at  sea  in  the  RBTS  mode,  given  a  sufficient  sam¬ 
ple  size. 

For  the  particular  COTS  configuration  studied,  it  was  determined 
that  any  requirement  for  additional  offload  capability  should  be  ob¬ 
tained  through  the  use  of  motion  compensating  cranes.  Offload  capa¬ 
bility  that  could  not  be  provided  through  the  utilization  of  the  max¬ 
imum  motion  compensation  capability  could  be  obtained  by  the  addition 
of  breakwaters  as  required.  Two  basic  COTS  designs  emerged  as  poten¬ 
tially  viable  candidates  from  the  optimization  of  cost  effectiveness: 
These  were: 

a)  Maximum  motion  compensation  gain  and  no  breakwater,  and 

b)  Maximum  motion  compensation  gain  and  a  breakwater  with 

a  cutoff  frequency  near  the  frequency  of  the  peak  energy 
density  of  the  sea  spectrum. 

5.2  RECOMMENDATIONS  FOR  FURTHER  RESEARCH 

This  methodology  relies  heavily  upon  Gaussian  inputs,  linearizable 
dynamic  properties  of  the  system  elements,  and  a  minimum  of  coupling 
between  the  dynamic  elements.  Future  research  should  examine  the 
possibility  of  handling  non-Gaussian  inputs,  system  dynamics  that  are 
not  readily  linearized,  and  systems  where  strong  coupling  exists  bet¬ 
ween  the  dynamic  elements. 
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In  the  particular  case  of  the  COTS  problem,  several  areas  warrant 
further  investigation: 

a)  More  complex  forms  of  motion  compensating  crane  strategies 
should  be  examined,  such  as  those  that  track  the  relative 
displacement  between  the  cargo  and  the  adjacent  craft  or 
optimal  strategies  that  anticipate  relative  motion  and 
land  the  cargo  on  the  peaks  where  relative  velocity  is 
minimized. 

b)  The  interactions  between  the  floating  elements  should  be 
examined  in  greater  detail.  Particular  attention  should 
be  given  to  the  sheltering  of  one  craft  by  another  and 
the  effects  of  secondary  wave  generation  from  the  motions 
of  the  larger  vessels. 

c)  Finally,  additional  data  needs  to  be  collected  in  order  to 
more  accurately  establish  the  component  cost  models  and  to 
verify  or  develop  improved  expressions  for  the  container 
offload  cycle  times. 

Based  upon  the  predicted  performance  of  the  simplified  motion  com¬ 
pensation  system  examined  in  this  study,  it  is  recommended  that  a  proto¬ 
type  motion  compensation  system  using  this  platform  stabilization  strategy 
be  designed  and  tested  on  a  floating  crane.  This  will  provide  an  interim 
motion  compensation  capability  until  the  feasibility  of  the  more  complex 
systems  is  demonstrated. 
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APPENDIX  A.l 


STANDARD  CRANE  COST 

The  large  number  of  crawler  cranes  in  the  larger  capacities 
currently  in  service  in  the  commercial  construction  industry  make  them 
viable  candidates  for  use  in  the  COTS  application.  As  discussed  in 
Sections  2.2.4  and  3.2.4,  the  parameters  that  most  strongly  influence 
the  cost  of  a  candidate  crane  are  type,  lift  capacity,  and  reach  require¬ 
ments.  For  a  crawler  crane,  the  weight  of  the  lifted  load,  W,  herein 
after  referred  to  as  the  lift,  and  the  reach,  R,  can  be  used  to  approxi¬ 
mate  the  cost.  The  load  moment  will  be  defined  as  the  lift  times  the 
reach,  UR.  Commercially  available  crawler  cranes  have  load  moments  up 
to  about  7  x  106  lb-ft.  If  crawler  cranes  were  made  with  lift-reach 
capabilities  above  this,  they  would  most  likely  be  more  expensive  than 
specially  modified  cranes  like  the  Manitowoc  4100W  Ringer.  This  is  a 
standard  crawler  crane,  outfitted  with  an  extra  long  mast  and  heavy  duty 
boom,  mounted  in  the  center  of  a  36  foot  diameter  steel  ring.  Extra 
counterweights  are  supported  by  rollers  on  this  ring  behind  the  crane 
and  the  boom  base  is  supported  in  a  similar  fashion  in  front  of  the 
crane.  The  cost  of  this  type  of  crane,  installed,  on  a  ship,  is  approx¬ 
imately  1.5  million  dollars  and  it  has  lift-reach  capabilities  in  excess 
of  16  x  106  lb-ft. 

For  lift- reach  requirements  below  7  x  106  lb-ft,  where  standard 
crawler  cranes  provide  the  more  economical  alternative,  a  reasonable 
cost  model  is  required. 
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The  lift  capability  of  crawler  cranes  drops  off  sharply  with  in¬ 
creasing  reach.  Because  overturning  of  the  crane  is  one  of  the  primary 
rating  criteria  for  these  types  of  cranes,  the  plot  of  load  moment,  WR, 
versus  reach,  R,  should  be  close  to  a  constant.  A  plot  of  load  moment 
versus  reach  for  three  different  commercially  available  cranes  is  shown 
in  Figure  A.l.  They  slope  downward  because  at  longer  reaches,  the 
weight  of  the  boom  and  rigging  reduce  the  useful  capacity  of  the  crane. 
They  are  not  straight  lines  because  the  tipping  moments  are  not  counter¬ 
acted  at  the  center  of  rotation  but  at  some  point  between  this  and  the 
tips  of  the  crawlers.  Straight  line  approximations  are  useful,  however, 
for  reach  requirements  of  container  handling  operations  (from  50  to 
130  feet).  It  can  be  observed  that  the  larger  capacity  cranes  have 
more  negative  slopes  of  their  load  moment  curves.  If  these  slopes  are 
plotted  against  an  arbitrary  measure  of  capacity,  such  as  load  moment 
at  70  feet  reach,  M7Q,  Figure  A. 2,  there  appears  to  be  a  reasonable 
linear  relationship. 

S  =  -.00463  M70, 

where  S  is  the  slope  of  the  load  moment  curve. 

The  cost  of  cranes  varies  with  the  options  available  on  a  partic¬ 
ular  model,  the  economic  conditions  at  the  time  of  purchase  and  the  com¬ 
mission  structure  of  the  particular  dealer.  The  prices  of  several 
crawler  cranes  obtained  from  telephone  interviews  with  dealers  is 
plotted  with  a  ten  percent  confidence  interval  against  their  respective 
load  moments  at  seventy  feet,  M7Q  (Figure  A. 3).  A  reasonably  linear 
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relationship  is  found. 

COST  =  .18  M?q. 

The  linearized  load  moment  equations  can  be  written, 

WR  =  SR  +  B, 

where  the  "y  intercept",  B,  is  a  function  of  crane  capacity. 

It  follows  that, 

M?0  =  (-.00463  M70)  70  +  B. 

B  =  1 .3741  M7q. 

Therefore,  WR  =  (-.00463  M7Q)  R  +  1.3241  M7Q. 

This  can  be  solved  for  M7Q  as  a  function  of  lift  and  reach. 

WR 

M70  =  1.3241  -  .00463  R  * 

Since,  Cost  =  .18  M7q, 

Cost  =  .18  WR  =  $1000  WR 

T7S&]  -  ToSSSn?  7356  -  25.76  R  ' 

Crawler  crane  cost  is  now  defined  as  a  function  of  lift  and  reach. 


If  lift  and  reach  define  a  cost  greater  than  $1.5  million,  it  will  be 
assumed  that  the  ringer  crane  described  earlier  would  be  chosen,  other¬ 
wise  a  standard  crawler  crane  will  be  assumed. 


APPENDIX  A. 2 


MOTION  COMPENSATING  CRANES 

A  number  of  researchers  have  addressed  the  problems  associated 

with  the  development  of  cranes  that  will  automatically  compensate  for 

( 8 

relative  motion  between  two  vessels  when  transferring  loads  at  sea.  ’ 

9,10,11  ,38,39,40,41  ,42)  _.  ,  •  .  *  ..  .  ..... 

The  COTS  design  goal  of  operational  capability 

in  a  sea  state  3  results  in  significant  maximum  relative  velocities  in 
the  neighborhood  of  10  ft/sec  (3  m/sec).  For  container  weights  up 
to  35  tons,  peak  synchronization  horsepowers  above  1200  HP  will  be  re¬ 
quired. 

One  of  the  major  hurdles  to  be  overcome  is  the  design  of  a  motion 
compensating  crane  is  the  stabilization  of  the  feedback  control  system. 
The  control  loop  includes  the  dynamics  of  the  control  strategy,  the 
dynamics  of  the  winch  or  other  motion  compensating  machinery  (i.e., 
hydraulic  cylinders,  counterweights,  etc),  and  the  lightly  damped  struc- 
ural  resonance  of  tne  crane  itself. 

One  study  has  shown  that  this  problem  alone  is  a  formidable 
one  because  of  the  proximity  between  the  natural  structural  resonance 
of  the  crane  and  the  frequency  range  of  the  lighter  motions. 

The  design  requirements  of  a  motion  compensating  crane  include 
rapid  response  to  control  inputs,  machinery  dynamics  that  are  well  above 
the  bandwidth  of  the  expected  disturbances  and  relatively  insensitive 

to  load,  and  a  capability  of  transmitting  large  amounts  of  horsepower. 
These  requirements  are  best  met  by  hydrostatic  transmissions. 
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A  hydrostatic  transmission  consists  of  a  variable  displacement 
pump  driven  by  a  prime  mover,  directly  connected  by  a  closed  hydraulic 
loop  to  a  fixed  displacement  hydraulic  motor.  The  stroke  of  the  pistons 
in  the  pump  are  controlled  by  the  angle  of  a  swashplate,  which  can  be 
electrically  actuated.  This  lends  itself  readily  to  series  compensation 
stabilization  techniques.  A  typical  transfer  function  for  such  a  system 
is : 


(1  +  Ts)  (1  +  + 


where  x  =  output  speed 
u  *  control  input 
K  =  system  gain 

t  =  swashplate  motion  time  constant 
“n^  =  hydrualic  system  resonance 
C  *  damping  ratio  . 


The  term, 

(1  +  ts) , 

represents  the  lag  response  of  the  swashplate  position  and  the  term. 


o  ♦£*♦£>, 


represents  the  second  order  effects  due  to  the  masses  of  the  moving 
machinery  elements,  the  damping  characteristics  of  the  hydraulic  oil  and 
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the  compressibility  of  the  oil  trapped  between  the  pump  and  the  motor. 
Typical  values  for  these  constants  are: 

t  =  .3  sec 

uj  *  10  rad/sec 
nn 

5  3  .5 

The  stability  problem  can  be  more  readily  appreciated  when  a  sys¬ 
tem  block  diagram  is  examined. 


CONTROL 

KG(s) 

1 

1 

2c 

(1  ♦  rs)  0  ♦  ) 

1 

,  s* 

2c 

STRATEGY 

s 

'Vv  " 

c 

=  vertical  displacement  of  the  lighter 
Zc  s  vertical  displacement  of  the  cargo 

Zt  =  vertical  motion  of  crane  boom  tip  if  crane  were  a  rigid  body 

wn  =  natural  resonance  of  crane  structure 

Standard  feedback  control  strategies  such  as  proportional  or  pro¬ 
portional  plus  rate  for  the  above  system  can  only  be  stabilized  at  very 
low  loop  gains  or  by  a  complex  network  of  compensation  lags  and  leads 
resulting  in  conditional  stability.  Since  the  structural  resonance 
varies  with  both  the  load  and  boom  position,  a  computer  will  probably 
be  required  to  continually  fine  tune  the  system  in  order  to  assure 
global  stability. 


A  great  deal  more  research  must  be  conducted  before  the  total , 
relative  motion  compensation/adapti ve  touchdown  strategies,  currently 
under  investigation  are  developed  to  a  point  that  they  can  be  effect¬ 
ively  included  and  analyzed  in  this  present  tradeoff  study. 

Crane  platform  roll  is  the  greatest  contributer  to  the  total  dis¬ 
placement,  but,  because  of  its  relatively  low  frequency,  is  a  minor 
contributer  to  the  total  relative  velocity.  This  indicates  that  a 
motion  compensating  crane  that  works  to  stabilize  the  platform  may  be 
an  approach  requiring  significantly  less  horsepower.  A  simplified 
motion  compensating  crane  design  is  proposed  here  as  a  technologically 
feasible  alternative.  The  key  elements  of  the  proposed  scheme  are: 

a.  Compensate  for  crane  platform  motions  only, 

b.  Use  velocity  feedback  instead  of  position  control, 

c.  Do  not  close  the  loop  around  the  structural  resonance  of 
the  crane, 

d.  Take  maximum  advantage  of  the  adaptivity  and  anticipatory 
characteristics  of  the  human  operator. 

A  block  diagram  of  this  platform  stabilized  motion  compensating 
crane  is  shown  below. 


m 


Zg  =  motion  compensated,  effective  boom  tip  motion 

Zr  =  vertical  relative  motion  between  cargo  and  lighter. 

The  first  control  loop  reduces  the  effects  of  the  absolute  vertical 
motions  of  the  boom  tip.  If  the  feedback  in  this  loop  had  sufficient 
bandwidth  and  gain,  it  could  hold  the  cargo  still  in  inertial  space  as 
the  boom  tip  rose  and  fell,  leaving  the  crane  operator  free  to  focus 
his  attention  on  the  motion  of  the  lighter  only. 

The  transfer  function, 


represents  all  of  the  characteristics  of  the  human  operator  from  pre¬ 
diction  to  automatic  gain  control.  The  human  operator  closes  the  con¬ 
trol  loop  around  the  structural  resonance  of  the  crane.  This  is  nor¬ 
mally  done  in  conventional  crane  operation.  If  a  lighter  had  a  heave 
RAO  of  1 ,  it  would  have  a  significant  amplitude  somewhat  less  than 
3  feet  and  a  predominant  period  around  6  seconds,  in  sea  state  3.  It 
is  not  unrealistic  to  expect  the  human  operator  to  be  able  to  anticipate 
peaks  in  this  motion  and  land  the  cargo  accordingly. 

One  problem  in  implementing  this  strategy  is  obtaining  the  Zt 
signal  free  from  contamination  by  the  structural  resonance  of  the  crane. 
This  can  be  best  accomplished  by  measuring  the  motion  of  the  platform, 
assuming  the  crane  is  a  rigid  body  and  calculating  the  theoretical 
motion  of  the  boom  tip.  A  second  problem  is  stability.  Even  this  rel¬ 
atively  simple  system,  in  the  absence  of  compensation,  G(s)  is  unstable 
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for  significant  values  of  loop  gain,  K.  High  loop  gain  is  important  in 
this  system  because  the  low  frequency  motion  attenuation  is  approxi¬ 
mately  equal  to  1  . 

K+T 

A  simple  analysis  using  Bode  plots,  indicates  that  acceptable 
stability  can  be  achieved  with  gains  as  high  as  10  if  a  lag-lead,  lead- 
lag  series  compensation,  G(s),  is  chosen. 

(1  +  .3s)  (1  +  .Is) 

G(s)  =  (1  +  s)  (1  +  .01s) 

The  resulting  closed  loop  transfer  function  between  rigid  body 
boom  tip  motion  and  the  effective  boom  tip  motion  is, 

h  *  (1  +  s)  (1  ♦  .Is  +  .01s2)  (1  +  ,01s) _ 

Zt  (1  +  s)  (1  +  .Is  +  .01s2)  (1  +  .01s)  +  K  (1  +  .Is)’ 

The  hoisting  velocity,  is  the  difference  between  the  boom  tip 
velocity  and  the  effective  boom  tip  velocity. 

v  *  -•*.  • 

The  transfer  function  between  the  boom  tip  motion  and  the  hoisting 
velocity,  V,  is, 

V_  =  _ s  K  (1  +  .Is) _  _ 

Zt  (1  +  s)  (1  +  .Is  +  .01s2)  (1  +  .01s)  +  K  (1  +  .Is)’ 

The  significant  hoisting  velocity  and  the  weight  of  the  cargo  can 
be  used  to  determine  the  approximate  value  of  the  horsepower  required  to 
Implement  this  motion  compensating  strategy.  The  transfer  functions 
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given  in  this  Appendix  will  be  used  in  the  COTS  model  to  calculate  the 
motion  attenuation  characteristics  and  horsepower  requirements  for  this 
design  of  motion  compensating  cranes. 

A  plot  of  commercially  available  diesel  engine  costs  versus  horse¬ 
power  is  given  in  Figure  A. 4..  In  general,  the  size  and  cost  of  the 
major  power  transmission  components  of  the  motion  compensating  system 
tend  to  increase  with  horsepower  requirements.  The  cost  of  the  elect¬ 
ronic  controls,  sensing  elements  and  installation,  however,  will  tend 
to  remain  a  constant  for  a  range  of  horsepowers  as  these  items  do  not 
change  significantly. 

Based  on  the  procurement  costs  of  the  hydraulic  power  pack  and 
winch  system  for  the  prototype  RBTS,  a  cost  of  approximately  $400.00  per 
horsepower  was  chosen  for  the  variable  cost  of  the  motion  compensating 
hardware. 

Based  on  current  projected  motion  compensating  crane  cost  esti- 
(21 

mates,  '  a  fixed  cost  of  approximately  $250,000.00  appears  to  be  rea¬ 
sonable  at  this  time. 

This  leads  to  a  motion  compensation  cost  model  of  the  form, 

COST  =  $250,000  +  $400  HPs , 

where  HPs  is  the  significant  horsepower  as  defined  in  Section  3.2.4. 


APPENDIX  B 
PRODUCTIVITY 


B.l.  Minimum  Time  to  Position  a  Pendulum 
B.2.  Cycle  Time 
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APPENDIX  B.l 

MINIMUM  TIME  TO  POSITION  A  PENDULUM 


In  the  development  of  the  lost  time  expressions  in  the  productivity 
model  (Section  3.3.2),  an  analogy  is  made  between  the  crane  operator's 
task  of  positioning  the  cargo  and  the  idealized  problem  of  positioning 
a  simple  pendulum.  It  is  assumed  that  an  operator  might  perform  posi¬ 
tioning  tasks  in  a  suboptimal  manner  in  times  that  are  longer  than,  but 
approximately  proportional  to  the  minimum  time  obtained  by  the  optimal 
solution.  This  minimization  problem  forms  the  basis  for  the  use  of  vr 
in  the  lost  time  expressions  of  the  productivity  model. 

Given  an  ideal  pendulum  at  rest,  that  can  be  moved  in  a  horizontal 
line  at  its  fulcrum;  find  the  minimum  time  and  optimal  trajectory  that 
will  move  the  pendulus  mass  to  an  adjacent  point  at  a  state  of  rest  in 
the  minimum  time  if  the  maximum  velocity  of  the  fulcrum  is  limited. 


x(0)  =  x(0)  =  x(0)  =  0 
y(0)  =  0 

x(tf)  =  y(tf )  *  xf 
x(tf)  *  x(tf)  =  0 
|y(t)|  <  v 


x  =  displacement  of  pendulus  mass 
y  *  displacement  of  fulcrum 
l  *  pendulum  length 
M  »  pendulus  mass 
g  *  acceleration  of  gravity 
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Equation  of  Motion 

Mx  =  ^p-(y-x) 

x  +  f*  =  |y 


Let 


=  a 

l 


x  +  ui2x  =  aj2y 

Because  the  constraint  is  on  y  rather  than  y,  a  more  useful  form 
of  the  equation  is: 

x  +  u>zx  =  o)2y 


Find  y*(t)  that  minimizes 


State  Equations 

Let  x  *  xlf  x  =  x2»  x  *  x3 ,  y  =  u 

Xi  *  x2 
x2  =  x3 

X3  =  -U)2X2  +  u)2U 


Therefore 


1 


Hamiltonian 

H  =  1  +  PjX2  +  P2X3  “  P3W^X2  +  P3U)^U 

Co-State  Equations 

f>i  =  0 

p2  =  -P1  +  u)2P3 


P3  =  P2 


Solution  of  the  Co-State  Equations 
P3  »  constant  *  Kj 

P3  =  -p2  =  P1  -  u>2p3 

P3  +  U2P3  =  Ki 

p3  *  k2  sin  (ait  +  $)  +  — 

U)2 

where  Klt  K2,  and  *  are  constants  of  integration. 
Optimal  Control 

u*(t)  is  the  value  of  u(t)  that  minimizes  the  Hamiltonian. 
Therefore, 
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u*(t)  =  -sgn{P3)  V. 

j f  %  <  M  ,  the  optimal  control  is  of  the  form 

T  —  U) 

(+V,  -V,  +V}  . 

Solution  of  Differential  Equation 

7  +  (u2x  =  +  tu2V 

x(t)  =  A  sin(wt)  +  B  cos  (wt)  +  C  +  Vt 

Let  tx,  t2,  and  t3  be  the  length  of  time  that  the  +V,  -V,  and  +V 
control  is  applied  respectively.  Substituting  these  times  into  the 
above  solution  to  the  differential  equation,  using  the  initial  and 
final  conditions  of  the  original  problem,  and  using  final  conditions 
of  each  phase  for  the  initial  conditions  of  the  next,  the  following 
set  of  algebraic  equations  result. 

sin  ( cot !  +  ut2  +  wt 3)  -  2  sin  (uit2  +  wt3)  +  2  sin  (iot3)  *  0 

cos  (ut3  +  wt2  +  <ot3)  -  2  cos  (wt2  +  wt3)  +  2  cos  (u>t3)  *  1 
xf 

ti  -  t2  +  t3  *  -f 

The  following  substitution  of  variables, 

Z3  =  ait  1  +  U)t2  +  o»t3 

Z2  =  (ot2  +  wt3 


z3  =  wt3 


results  in  a  set  of  equations  which  lend  themselves  nicely  to  Newton- 
Raphson  techniques  of  solution. 

sin  ( z ! )  -  2  sin  (z2)  +  2  sin  (z3)  =  0  . 

cos  ( z j )  -  2  cos  (z2)  +  2  cos  (z3)  =  1  . 

WX  » 

zx  -  2  z2  +  2  z3  =  -y- 

Note  that 


Appl i cation 

The  above  technique  was  used  to  solve  for  the  minimum  time  to  re¬ 
position  a  pendulum  a  distance  of  two  feet  at  a  maximum  speed  of  two 
feet/second  for  a  range  of  pendulum  lengths.  The  basic  computer  pro¬ 
gram  used  to  implement  this  analysis  is  given  in  Figure  B.l.  The  re¬ 
sults  are  plotted  in  Figure  B.2. 

Values  of  pendulum  length  between  14  feet  and  130  feet  are  plotted. 
It  can  be  seen  that, 

T  =  1  +  .  36  , 

provides  a  reasonable  fit  to  these  data.  Since  the  minimum  time  to  re¬ 
position  a  rigid  body  two  feet  at  a  maximum  speed  of  two  feet/second  is 
one  second,  the  lost  time  attributed  to  the  pendulum  is 

Tlost  ■  -36^- 

In  general,  therefore,  for  pendulum  movements  that  are  small  with 
respect  to  the  pendulum  length,  it  can  be  shown  that  the  extra  time 
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Figure  B.l.  Basic  Program  Used  to  Find  the  Minimum  Time  to 
Reposition  a  Pendulum  a  Distance  of  Two  Feet  at  a  Speed 
of  2  Feet/Second  for  a  Range  of  Pendulum  Lengths 
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(lost  time)  required  to  complete  a  positioning  task  is  approximately 
proportional  to  the  square  root  of  pendulum  length. 


APPENDIX  B.2 


CYCLE  TIME 


An  estimate  of  the  container  offload  cycle  time  is  required  as 
part  of  the  productivity  model  developed  in  Section  3.3.2.  A  cycle 
time  model  of  the  form. 


CT 


T_  +  T  +  L 


» 


was  chosen,  where 

CT  =  cycle  time, 

T  =  minimum  time  to  move  a  container 
m 

under  ideal  conditions, 

T$  *  average  time  lost  over  the  ship, 

T^  =  average  time  lost  over  the  lighter. 


It  was  hypothesized  that  the  lost  times  could  be  approximated  by  a 
power  function  of  the  two  parameters  that  most  directly  influenced  it  - 
relative  vertical  displacement  between  the  cargo  and  the  adjacent  vessel 
and  pendulum  length  beneath  the  crane. 

T  *  A 

where 

T  =  lost  time 
A  =  constant 

o  *  standard  deviation  of  vertical  displacement 
l  =  pendulum  length 
a  and  8  *  exponents 
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After  the  OSDOC  II  exercises,  several  observers,  knowledgeable  in 
the  field  of  marine  operations  and  Naval  material  handling,  extrapo¬ 
lated  the  data  to  generate  estimates  of  offload  cycle  times  for  sea 
states  1  through  4.^)  j^g  minimum  time,  T  was  subtracted  from 

these  estimates  to  obtain  the  total  lost  time.  A  review  of  the  data 
indicated  that  the  time  lost  over  the  ship,  T  ,  was  approximately  70% 
greater  than  the  time  lost  over  the  lighter,  T^.^  This  data  is  pre¬ 
sented  in  Table  B.l.  along  with  estimates  of  the  standard  deviation  of 
relative  vertical  displacements  over  the  ship  and  lighter  for  the  re¬ 
spective  sea  states.  The  vertical  displacements  were  obtained  from 
analysis  of  a  typical  container  offloading  configuration  by  the  RELMO 
computer  program. 

These  cycle  times  were  for  standard  cranes,  and  therefore,  the 
pendulum  lengths  were  assumed  to  be  120  feet  over  the  ship  and  150  feet 
over  the  lighter.  In  Appendix  B.l,  it  was  shown  that  the  exponent,  /3, 
could  be  taken  as  h- 

There  are  now  only  two  unknowns  in  each  equation  -  A  and  a.  A 
log-log  plot  of  the  data  in  Table  B.l  (Figures  B.3  and  B.4)  allows 
these  remaining  two  unknowns  to  be  approximated. 

The  resulting  lost  time  expressions  over  the  ship  and  the  lighter 
respectively  are: 

Ts  ■  ■8iVs‘s 
h  ■  tli  ■ 


and 


TABLE  B.l.  DATA  FOR  EVALUATION  OF  SHIP  AND  LIGHTER  LOST  TIME  EXPRESSIONS 
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(4)  Results  of  application  of  RELMO  Computer  program  to  typical  configurations. 
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APPENDIX  C 


TETHERED  FLOAT  BREAKWATER 

This  Appendix  contains  correspondence  between  the  author  and 
faculty  at  Scripps  Institution  of  Oceanography  regarding  cost  and  per¬ 
formance  of  the  tethered  float  breakwater.  The  data  contained  herein 
were  used  to  develop  the  functional  relationships  given  in  Section 
3.2.2  and  used  in  the  TFB  subroutine  of  the  COTS  computer  model  dis¬ 
cussed  in  Section  3.4. 
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CodeAQZ2 
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„  STATE  Of  CALIFORNIA 


UNITED  STATES  ARMY 
Corps  of  Engneers 
MAR/T1ME  ADMMSTRATION 

February  25,  1977 


Mr.  Dextar  Bird 

•-.'ashing ton  Analytical  Services 
Cantar,  Inc. 

IG&G 

1130  Fields  Road 
Aockvllle ,  Maryland  20850 

Dear  Daxtar: 

It  was  a  pleasure  to  meet  with  you  and  Las  this  week  and 
to  discuss  the  Crana/TFB  Optimization  Study.  I  have 
given  it  some  considerable  thought  since  you  left  and 
1  feel  that  Che  approach  you  are  talcing  is  quite  sound 
and  should  result  in  a  good  feel  for  the  optimum  system. 

When  you  have  decided  on  your  example  sea  states  (or 
representative  Pierson-Moskowitz  spectra),  send  them 
along  and  I  will  respond  with  families  of  transfer 
functions  as  ve  discussed  and  with  costing  formulas. 

3est  of  luck  on  your  work. 

Sincerely, 


R.  J.  Seymour 


cc:  Mr.  Leslie  Bonds 

Mr.  M.  £.  Essoglou 


I 
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I 

I 
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J^EGcG 


WASHINGTON  ANALYTICAL  SCNVlCCS  CINTf  A.  INC. 

H  V  DRO  SPAC  f  *CHALLC  NGC  R  GJtOU* 

2150  P'JArt 
RockvHle.  Maryland  20650 
(301 )  W8-4J50 

J-5772-1 5 
29  March  1977 
JW8:nbt 


Or.  R.  J.  Seymour 

The  Tethered  Float  Breakwater  Ocean  Experiment 
Mail  Code  A033 

ucso 

La  Jolla,  California  92093 


Re:  Breakwater  Characteristics  and  Costs 
Dear  Sir: 

In  accordance  with  your  recommendations,  I  am  beginning  the 
Brea kwater/Mot Ion  Compensation  trade-off  analysis  based  upon  sea  state 
3  only. 

It  is  most  convenient  for  me  to  use  the  following  formulation  of 
the  Plerson-Moskowltz  sea  spectrum: 


S(  ui  )  »  8.4  exp 

T 

Hi 

where  *  significant  wave  height  in  ft. 

»  *  angular  frequency  In  rad/sec 
and  S  (  at  )  duj  * 

2 

where  r  »  variance  of  sea  level  distribution  about  the  mean 

In  ft2. 

1  would  like  to  define  sea  state  3  as  having  a  significant  wave  height 
of  5  ft.  -This  yields  a  sea  spectrum  with  a  variance  of  approximately 
1.5625  ft\ 


-33.6 

HTT 

uj 
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29  March  1977 
Page  2 


As  we  discussed  earlier,  I  need  some  information  on  the  Energy 
Transmission  Ratio  (ETR)  and  the  related  breakwater  costs.  For  ex¬ 
ample,  if  the  ETR  can  be  described  by  a  function  of  the  form: 


ETR  »  ETR  (a»c,  a  ), 


it  would  be  preferable  to  have  a  cost  function  of  the  form: 

Cost  *  C8  (»c,  a ). 

If  this  formulation  introduces  unrealistic  constraints,  please  let 
me  know  and  we  will  try  another  approach. 

Thank  you  again  for  your  continuing  assistance  in  this  effort. 


Cordially, 

EG4G  WASHINGTON  ANALYTICAL 
SERVICES  CENTER,  INC. 


Copy  to: 

Mr.  M.  E.  Essoglou,  NAVFAC  031 A 
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May  20,  1977 


itr.  Oexter  9ird 

Jashington  Analyti  :at  Services 
Center,  Inc. 

2G&G 

2150  Fields  Road 
Rockville,  Maryland  20850 

Dear  Mr.  Bird: 

The  performance  and  cost  of  a  tethered  float  breakwater  for 
a  given  sea  state  are  related  to  several  design  parameters. 

For  example,  float  characteristics,  tether  length,  and  the 
number  pf  rows  are  all  significant.  However,  all  functional 
relationships  which  exist  between  these  variables  are  highly 
nonlinear,  which  makes  even  a  simple  analysis  somewhat 
compl icated. 

I  have  investigated  the  performance  and  costs  for  several 
different  breakwater  designs  subject  to  sea  state  three, 
with  a  significant  wave  height  of  5.2  feet,  in  water  of 
depth  equal  to  50  feet.  For  each  design  I  have  enclosed 
a  plot  of  energy  transmission  ratio  (ETR)  vs.  Frequency  (j), 
and  an  estimate  of  the  cost  per  foot  of  breakwater  provided 
by  Oick  Seymour.  ETR(w)  is  defined  as  the  ratio  of  trans¬ 
mitted  energy  at  frequency  w  to  incident  energy  at  frequency  w. 

There  are  several  bits  of  information  contained  on  each  plot. 
They  have  the  following  meanings:  Teth  is  the  effective 
tether  length  in  cm,  o  is  the  specific  gravity  of  the  float, 

'1  is  the  number  of  rows  of  breakwater,  and  ETR  overall  is  the 
ratio  of  total  energy  transmitted  to  total  energy  incident. 
Please  ignore  the  line  labled  "measured. “  Keep  in  mind  that 
a  value  of  ETR  at  a  certain  frequency  has  meaning  if  there  is 
a  reasonable  amount  of  incident  energy  at  that  frequency. 

For  the  sea  state  considered,  the  range  of  incident  energy 
is  approx imatel y  from  0.08  to  0.5  hertz. 
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For  some  sea  states  and  breakwater  designs,  the  TF8  functions 
very  mucn  like  a  low  pass  filter.  For  such  a  system,  per¬ 
formance  can  be  approximated  by  a  function  of  the  form: 


*nere: 

(1)  ETR  *  1  +  .  95/m)H(w-a,c- .  95/m)  -  m(w-«c  )H(w*ac  )  . 


H(arg) 


s(0  for  arg<0 
(1  for  arg>0 


dc  and  m  are  defined  in  diagram. 


It  appears  that  this 
des i gns : 

approximation 

is  valid 

(  -1> 
m(hz  ) 

for  the  following 

Cost 

Design  ♦ 

-c(hz) 

(S/Foot) 

1 

.306 

-1  3 

300 

3 

.275 

-13 

750 

S 

.244 

-1  3 

1050 

10 

.213 

-14 

2250 

14 

.205 

-13 

2400 

1  7 

.197 

-13 

3000 

21 

.1  90 

-1  3 

3000 

25 

.182 

-13 

3000 

29 

.  1  60 

-1  6 

5250 
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I  nave  attempted  to  express  cost  as  a  function  of  j 
I  found 

cost  :  5.4  j,.*3'8 

$/foot  j  in  (hz) 

for  .  306^c^_.  150  fiz 

and  for  »  :  0.05 

The  other  designs  (with  a>0.05)  produce  curves  which  do  not 
lend  themselves  to  simple  approximations.  So  I  have  included 
copies  of  all  of  the  designs  with  their  associated  costs. 

You  may  be  aole  to  construct  functions  for  these  that  will 
meet  your  needs. 

let  me  know  if  there  is  any  additional  information  I  can 
provide, (714)  452-2561  . 

Very  truly  yours, 


Daniel  Hanes 


Enclosures 

cc:  Mr.  Don  B.  Jones 


UNIVERSITY  OF  CALIFORNIA.  SAN  DIEGO 


UNIVERSITY  OF  CALIFORNIA.  SAN  DIEGO 
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APPENDIX  D 


EVALUATION  OF  RAPS 

Program  Description 

This  appendix  describes  the  input  data  required  by  the  RAOS  program 
given  in  Appendix  E  in  order  to  generate  the  response  amplitude  oper¬ 
ators  of  the  floating  vessels.  The  RAOS  program  is  a  modified  version 
of  RELMO,  a  computer  program  developed  by  the  Navy  to  solve  for  relative 
motion  between  ships  in  random  seas.  A  complete  description  of  the 
original  RELMO  program  can  be  found  in  Technical  Note  N-1371,  The  Motion 
of  Floating  Advanced  Base  Components  in  Shoal  Water  —  A  Comparison 
Between  Theory  and  Field  Test  Data  by  D.A.  Davis  and  H.S.  Zwibel  of  the 
Civil  Engineering  Laboratory,  Port  Hueneme,  California. 

The  modified  program,  RAOS,  makes  use  of  that  portion  of  RELMO 
necessary  to  evaluate  the  vessels  response  amplitude  operators  and 
writes  them  on  disc  files  to  be  recalled  later  by  the  COTS  program  (also 
listed  in  Appendix  E).  The  RAOS  program  consists  of  a  MAIN  program  that 
formats  the  data  output  and  four  subroutines  as  follows: 

RAO:  Reads  input  data  and  calls  other  subroutines. 

FREQ:  Generates  angular  frequencies  and  wave  numbers  to  be  used 

In  computing  frequency-dependent  terms  In  equations  of 
motion.  Computes  associated  depth-dependent  wave  lengths. 

ADMAB:  Computes  added-mass  and  damping  coefficients  by  Grim's 
method . 
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COEF:  Computes  force  coefficients  appearing  in  equation  of 

motion.  Dynamic  force  coefficients  are  functions  of  added- 
mass  and  damping  coefficients  computed  in  subroutine  ADMAB. 
Solves  for  the  frequency-dependent  ship  coordinates  (heave, 
pitch,  surge,  yaw,  roll,  and  sway).  These  are  the  response 
amplitude  operators  for  each  ship  as  all  exciting  waves  are 
assumed  to  have  unit  amplitudes.  The  process  involves  in¬ 
version  of  the  force  and  moment  coefficient  matrix  and 
multiplication  of  inverted  matrix  with  the  wave  force  and 
moment  matrix. 

Program  Inputs 

Data  are  entered  on  two  types  of  input  cards:  INPUT  and  SHIP. 

INPUT:  Data  needed  to  define  the  period,  length,  and  direction  of 
the  unit  amplitude,  incident  regular  waves  is  placed  on 
this  card.  Card  variables  are  defined  as  follows: 

DEPTH:  The  mean  water  depth,  in  feet,  at  a  station  beneath 
the  center  of  gravity  of  the  ship. 

NFREQ:  The  number  of  angular  frequencies  to  be  used  in 
computing  the  ship  response  amplitude  operators. 

This  variable  and  the  one  designating  the  number  of 
ship  transverse  sections  (NSEC)  have  a  predominant 
effect  on  the  amount  of  computational  time  required. 
Computional  time  will  increase  roughly  as  the  pro¬ 
duct  of  NFREQ  and  NSEC.  NFREQ,  however,  should  be 
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at  least  50  to  insure  sufficient  accuracy  in  most 
ship  motion  problems. 

FREMIN:  The  minimum  angular  wave  frequency  (rad/sec) 
desired. 

FREMAX:  The  maximum  angular  wave  frequency  (rad/sec) 

desired.  Intermediate  wave  frequencies,  totalling 
NFREQ-2  in  number,  will  be  spaced  equally  between 
the  extremes  defined  by  FREMIN  and  FREMAX. 

BETA:  The  angle  of  approach  (degrees)  of  the  Incident 

unit  amplitude  waves.  BETA  equal  to  zero  de¬ 
signates  waves  approaching  sternon  (see  figure 
below). 


ISTOP:  This  is  a  control  variable.  For  any  ISTOP  value 

(integer)  other  than  zero,  the  computer  reads  a 
new  set  of  data  cards  (beginning  with  INPUT)  and 
repeats  an  entire  set  of  computations. 

SHIP:  Data  concerning  ship  characteristics  are  Inputed  on  this 

card.  All  weights  are  in  pounds,  and  distances  are  in 


feet. 
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ISHIP:  Used  by  RELMO  to  designate  ship  1  or  2.  In  the 

modified  version,  RAOS,  it  Is  always  assigned  a 
value  of  1 . 

NSEC:  The  number  of  transverse  ship  sections.  There 

can  be  as  few  as  one,  but  the  usual  number  is  10 
or  20. 

W:  The  weight  (displacement)  of  the  ship  (lbs). 

BG:  The  vertical  separation  between  the  ship's  center 

of  gravity  (CG)  and  center  of  buoyancy  (CB).  BG 
is  positive  if  CG  lies  above  CB. 

RG:  The  radius  of  gyration  of  the  ship  in  pitch.  If 

unknown,  it  should  be  set  equal  to  1/4  of  the  ship 
length. 

XCGFS:  The  logitudlnal  distance  of  the  ship  CG  from  the 

aft  perpendicular. 

B:  The  ship  beam  at  the  waterline  for  each  transverse 

secti on. 

SECOE:  The  ship  section  coefficient  at  each  transverse 

section.  Note:  Bulbous  sections  will  have  co¬ 
efficients  greater  than  unity.  Grim's  method  for 
computing  the  section  added-mass  and  damping  will 
provide  rather  poor  estimates  if  SECOE  greatly 
exceeds  one. 

DRAFT:  The  vertical  distance  from  the  waterline  to  the 

keel  at  each  transverse  section. 


SECL:  The  thickness  of  each  transverse  section.  Usually 

the  thinkness  is  the  same  for  all  sections. 

KX:  The  effective  surge  mooring  spring  constant  in 

units  of  Ib/ft,  assumed  to  be  0.0  if  not  entered 
on  SHIP  card. 

ZC8M:  The  vertical  center  of  buoyancy  at  each  transverse 

section. 

TPHI:  The  vessel  natural  roll  period  in  seconds.  If 

unknown,  the  roll  period  can  be  estimated  from 
the  following  expression: 

™ 

where  B  is  the  beam  of  the  vessel,  GMT  is  the 
transverse  metacentric  height,  and  C  is  an 
empirical  constant  which  varies  from  0.44  for 
full -formed  merchant  ships  to  0.39  for  more  round 
hull  ships. 

0G:  The  vertical  distance  between  the  ship's  CG  and 

the  water  plane.  OG  is  positive  if  the  CG  lies 
above  the  water  plane. 

RGZ:  The  radius  of  gyration  of  the  vessel  in  yaw.  If 

unknown,  it  should  be  set  equal  to  the  radius  of 
gyration  in  pitch. 


ZETA:  The  roll  damping  coefficient.  Use  ZETA  =0.04  for 

round  hull  vessels  and  ZETA  =  0.08  for  vessels 
with  flat  bottom  hulls. 

The  input  data  used  to  generate  the  RAO's  for  the  example  system 
configuration  described  in  Section  4.1  are  given  below: 


•  Heavily  loaded  C7  container  ship: 


♦  INPUT  DEPTH=*30.  *NFREG*80*FREHIN« .03*FREMAX-4.  *BETAD-133.  *ISTOP*0  • 

•SHIP  ISHIP»l*NSEC-20*W=6.502E+7* BG-16* 1 *RG-176. *XCGFS-326. *NOMORE=l * 

B(1 > -20. 3 » 40. 5 *55. 3 » 63.8.72.9.  78. 1 *82.1 * 84 . 4  *  85 . 0  *  05 . 0 , 83 . 0  *  04 . 0* 81 • 1 » 76.4* 
69. 7*60*2*47.1*31.. 6*16. 6*4 . 8  *  SECOE< 1 )«  »  394  * . 359  * . SOB  * .63V  * . 753  » . 838 * ,893, 
.937* .965* .977* .976* .968* .958* .941* . 906  * .856 . . 813  * . 792  * .843* 1 . 000* 

DRAFT (1 >*18.0*31.2* 30. 9 *30. 6 *30.4*30.1*29. 8 *29. 5 *29. 3* 29. 1*28. 8* 28. 3* 20. 2* 
27.9*27.7*27.4*27.2*27.0*26.8* 13.4, SECL<  1  >-20*33  .S.ZCBHU  > *9. 0.15. 6* 15. 5* 
15.3*15.2*15.1*14.9*14.8*14.7*14.6*14.4*14.3*14.1*14.0*13.9*13.7*13.6*13.5* 
13.4*13.2*  TPH I “30. 2* GMT* 1 . 33*0G*2. 9*RGZ-176. * ZETA-. 06  * 


•  Medium  loaded  C4  container  ship: 


•INPUT  r>EPrH=50. *NFREa=80*FREMIN».05*FRENAX«4, ,BETAD»135. *IST0P»0  • 

•SHIP  ISHIP*l*NSEC*20*W«3.624£+7.BGM2.4*RG=»I43.3*XCGFS-287.0*NOM0RE*0* 

B<1 >-5.6* 19. 8* 33. 8 *46. 5 *61. 7 *71. 7 *78. 2* 81. 3 *3*82. 0*81. 4 *77. 7* 70. 0*59. 4* 47.0* 
34.1*21.7*11*2*3.4*  SECOE ( 1 >  * . 400  * .453* .503* .553* .612* .689. .775* .055* .917* .946* 
.933* .893* .846* .808* . 779* . 750* . 753  * . 798* 1 . 0* 1 . 0* PRAFT< 1 > *  12 . 3* 19*24 . 6* 

SECL < 1 ) -20427 » 225  *  ZCBM ( 1 ) -6 .2* 19*12 • 3*  TPHI -12. 4  *GMT  *8 .0*00-7 • 4  *RGZ-1 43 . 3* 
ZETA-. 06  • 


•  LCM8,  landing  craft: 


•INPUT  DEPTH-50. *NFREO-0O*FR£HIN- ,05*FREMAX-4. *BETAD-133. * ISTOP-O  * 

•SHIP  ISHIP- 1  *  NSEC- 1 7 *  U- 1 . 15E+5 . BG-4 .13*  RG- 17.89*  XCGFS-34 . 3  *  NOHORE-O  * 

BC1 >-21 . *21. I* 12*21.2*21 .1*20. 9, 20.6* SECOE ( 1 >-l . *  1 . * . 982* . 953* . 920* .889* 
.066*4*. 850*. 847* . 840* .826* .005* . 767* . 681  *  DRAFT < 1 >- . 26* . 76* 1 .24  *  1 .66* l . 98, 
2.20*2.34*4*2.50*2.35*2.22*2.05*1.83*1 .51*1 .01 *SECL<1 >- 17*3.0* ZCBH(1)«. 13* 
.38* .62* .83* .99*1,10*1.17*4*1 .25* 1.18* 1*11 *1.03*. 92* . 76* . 31  * TPHI-3 .0* 
GMT-18.3  *00-3. 2*RGZ-17.09*  ZETA- • 08  *KX-0.  * 
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INSERT  I 

Changes  To  The  COTS  Program  Required  To 
Evaluate  Effect  of  Breakwaters  (Section  4.3) 
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00  90  I  *  :t.2 
MCO  -  I~1 
DO  80  J  ~  1.8 
J1  =  9— J 

BCF  =  FLOAT  (Jl)/2, 

CALL  EVAL  (TCCOST.  BUCOST.  T SCO ST ,  OLR ) 

WRITE  (A? 50)  BCF »  S2»  BWCOST*  TCCOSTf  TSCOSTe  OLR 
FORMAT  (IHj  6F12 ♦ 2/ > 

CONTINUE 

CONTINUE 


INSERT  II 

Changes  To  The  COTS  Program  Required  To 
Implement  The  Two  Parameter  Evaluation  Technique  (Section  4.4) 
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MCO-2 

NC-?. 

D  D  9  0  I  —  1  *  6 

FK-2 *  AFLOAT < 1 1 ) 

DO  30  J~l»10 
J 1 - 1 1  ~  J 

BCF- , 75+ . 05* FLOAT ( J1 ) 

B 1  =» .  2  5  *  F  L  0  A  T  <  J 1  >  -  ♦  2  5 
IF ( J1. . GT ♦  5 )  B'CF-Bl 

C  A  L  L  E  0  A  L  <  T  C  C  0  S  T  *  B  W  C  0  3  T ,  T  3  C  0  3  T  ?  0  L.  R ) 

Xl-FK 

X2-BCF 

X3-S2 

X4-BWC03T 

X5-TCC0ST 

X 6 -TSCOST 
X7-0LR 

X  8 -  T  3  COST  /  0  L.  R 
FK-l  ♦  l'*FK+  ♦  1 

CALL  E'v'AL.  ( TO  CO  ST  t  BWCOST  *  TSCOST » OUR ) 

X  9  —  (  TSCOST — X  6  >  /  ( 0  L  R  ~  X  7 ) 

FK--X.1 

BCF- . 9* BCF 

CALL  LVAL  (  1  CCOS  I  *  BWCOS  T  r  TSCOST  v  UI..R  ) 

X 1 0-  <  TSCOST-Xo )  /  <  0L.R-X7  > 

V  i  ■  .i.  r  E  6  *  4  9  )  X 1  ?  X  2  ?  X  ■■■•  >  a  4  ’>  x  ?  X  6  x  ,  '  >  3  » X  9  >  X 0 
F  OEM  A  T  ( 1H  ?  or  6  .>  2  ’  3F  3  ■■  0  4F3  *  2  ? /  > 


INSERT  III 

Changes  To  The  COTS  Program  Required  To 
Minimize  Cost  Effectiveness  (Section  4.4) 
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('! I  T  -  0 
riCO-2 
Nt;«2 

X  0  ~  1  4  15 

Y  0—0  ■>  0 
T  AU-  2 
EPS- *  02 
BCF— X0 
FK-YO 

CAL  L  EMAL < TCCOST > BWOQST , IS COST  * OL 
Z 1 ~  B  C  F' 

Z2:-FK 

Z3~T8C0ST/0LR 
i"  K  1  *  I.  X  F  K  f « 0  5 

C  A  I. .  L.  E  U  A  L  <  T  C  C  0  S  T » B  U  C  0  S  T  ,  T  S  C  0  S  T  » 0 1... 
Z4~ < TSC0ST/0LR-Z3 ) / < FK-Z2 > 

FK-Z2 

BCF-  4  98XBCF 

C  A  L.  L  l£  V  A  L  ( T  C  C  0  S  L  » B  W  C  0  S  T «  T  S  C  0  S  T » 0 1... 
Z5~ <  TSCOST /QLR-Z3 > / < BCF-Z 1 > 

Z  6  ■  S  Q  R  T  <  2. 4  %  Z  4  r  Z  5  :*  Z  5  > 

ST3-TAU 
NI  T-'NI  I  v- 1. 

IF  (NIT » GT,  1.000)  GO  TO  40 

X1=X0~STS*Z5/Z6 

Y 1  —  Y  0  •••■  3  f  S  #  Z  4  /  7. 6 

1. 1"  v  A 1  *  0  I  ♦  4  *  0  )  a  I  -4  *  v 

I F  ( X 1 «  LT  »  0  •>  0  )  X 1  -0  »  0 

.1.  F  ( Y 1  >  G  T  4 1 0 » 0 )  Y I  —  1 0  4  0 

IF  <  Y1  ♦ LT  *  0  4  0 )  Y 1-0.0 

A  S  T  S  -  S  Q  R  T  <  <  X 1  -  X  0  >  *  *  2  +  <  Y 1  -  Y  0  )  X  X  2 ) 

IF  (  A8TS  4  LT  4  EPS  '■  GO  TO  40 

BCF-X1 

FK-Yl 

Call  e val  (  t  ccost  *  bwcos  t  »  t  scost  »  ol 

Z7-TSC0ST/0LR 

IF ( Z7 4 GE 4 Z3)  GO  TO  30 

XO-Xl 
YO-Y  :l 
GO  TG  10 

3  i  3 — 3  i  3  /  2  * 

GO  TO  20 

!4P  .1.  r !::  %  ■"!  ■  GO  *  X 0  v  Y 0  v  Z.  3 "  N  i.  I 
OPH A  T  ■  1. H  4  3  >:  3  <.  >  F  i.  2 . 4  '•  •  .'■■■  ■'  i  .'.I  • 
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INSERT  IV 

Changes  To  The  COTS  Program  Required 
To  Implement  Constrained  Optimization  (Section  4.4) 


I 
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GLfiM- 1 70 
hCO  2  . 

NO  2 
X0=1 .5 

Y  0=1 0.0 
TAU- ,  2 
EF'S- ,  0002 
10  BCF-XO 
FK-YO 

CALL  EOAL < TCCOST . BWCOST ?  TSCOST ? OLR > 
Z 1  -  B  C  F 

22= FK 

Z  3  ~  I S  COST  f 5 0 0  ♦  *  ( 0  L  R  M  /  0 1...  R  >  *  *  5  0  0 
FK=1 . 1*FK+ *  05 

CALL  E V A L  <  T C C 0 ST  y  B W C 0 S T  r  T 3 C 0 3 T ,  0 1. R  ) 
'  ZOTSCOST+3 0 0 ,  *  ( 0 L  R  M  /  0  L  R )  *  *  S  0  0 

Z4«(Z4~Z3>/ ( FK-Z2 > 

FK=Z2 

BCF- ♦ 995*BCF 

CALL  EVAL t TCCOST » BWCOST - TSCOS T t OLR  > 
25=  ( TSC03T/0L. R - Z 3 )  /  ( B C F - Z  I.  ) 

Z 5= TSCOS Tf 500 . * ( OLRM/OLR > **500 
Z5= ( Z5-Z3 ) / <  BCF-Z 1  ) 

Z6-S 0 R  T ( Z  4 *  2  4  r Z  5  *  Z  5 ) 

s  rs= taij 

20  N  [  on  i:  T  -K 
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.[  F  (  N I T  »  0  r  ,  1 000  )  GO  TO  40 
X1--X0-3T  S  ♦  Z  5  /  7. 6 
Y 1 ~Y0~ 1 0 . #STS*Z4/Z6 
IF  (XI  .  G  'f  .4.0)  XI  =4.0 
I F ( X 1 .LT. 0 ♦ 0 )  X 1 -0 ♦ 0 

I F  <  Y 1 ♦ GT « 1 0 . 0  >  Y1-10.0 

IF  <  Y1  ♦ LT . 0 ♦ 0 )  Y 1-0 . 0 

ASTS=SGRT ( < X 1 -XO ) **2+ ( Y 1 - YO ) **2 > 

IF (ASTS.LT .EF'S)  GO  TO  40 

BCF-X1 

FK--Y1 

CALL  EVAL  <  TCCOST ?  BUCOST ? TSCOST . OLR ) 
Z7-TSC0ST  +  500 . *  <  OLRH/OL R ) Y^SOO 
IFCNIT.LT. 50)  GO  TO  25 

IFCZ7.GE.Z3)  GO  TO  30 
CONTINUE 

X0«X1 
Y0-Y1 
00  TO  10 
3TS-STS/2 . 

GO  TO  20 

UR I TE  C 6  9 50 )  XO  y  YO  >  Z3 1 N I T 
W R I T E ( 6  t 50 )  TSCOST 1 0LR y Z3 1 NI T 

BCF-XO 

FK=YO 

CALL  EVAL ( TCCOST  t BWGOST , TSCOST v OLR ) 
W  R I T  E  <6  y  5  0  )  TSCOST  y  0  L  R  *  Z  7  •>  N I T 
F  0  R  M  A  T  C 1 H  1 3  C  5  X  y  F  .1 2  ♦  4  )  y  8  X  y  1 5  ) 
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PROGRAM  RAPS 

Generates  Response  Amplitude  Operators 
For  The  Floating  Vessels 
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